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AN ANALYSIS OF A CHARRING ABLATION 


THERMAL PROTECTION SYSTEM 
By Donald M. Curry 

SUMMARY 


An analytical model is presented for predicting the transient one- 
dimensional thermal performance of a charring ablator heat protection 
system when exposed to a hyperthermal environment. The heat protection 
system is considered to consist of an ablation material and backup 
structure. The ablating material is further considered to consist of 
three distinct regions or zones ; char, reacting, and virgin material. 

A Fortran TV digital computer program (STAB II) utilizing an im- 
plicit finite difference formulation has been written for the IBM 7094/^0 
computer system. The program considers one ablating material and a 
maximum of 12 backup materials with conduction or radiation and/or con- 
vection allowed between materials. Thermal properties of all materials 
are temperature dependent with the properties of the charring material 
also state dependent. 

The governing differential equations and their implicit • finite dif- 
ference formulation are presented. The program input and output is de- 
scribed in detail. Also, a comparison of theoretical and experimental 
results is shown. 


INTRODUCTION 


The analysis and design of theimal protection systems for entry 
into an atmospheric environment has resulted in a voluminous amount of 
literature on the general subject of ablation. (See refs. 1 and 2 for 
a survey of information on ablation. ) The ablation materials may gen- 
erally be classified into three categories; subliming, melting and va- 
porizing, and charring. The charring ablator normally provides the 
most efficient thermal protection shield for the major portion of a 
manned entry vehicle. This report describes a method for predicting 
the thermal response of a typical charring ablation material. The 
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response of a charring material to a hyperthermal environment is ex- 
tremely complex and the mathematical model presented to analyze the 
transient behavior of the material contains simplifying assumptions and 
approxi m ations necessary to afford even a numerical solution. 

The equations derived in this analysis have been programed in 
Fortran IV for an IBM 7094/40 computer system. The numerical formula- 
tion of this digital program, designated STAB II, is such that an im- 
plicit solution is obtained. The input, output, and various program 
options are discussed in detail. 

The thermal response of a typical charring material as predicted 
by STAB II is compared with arc tunnel test results. As shown, the 
predicted in-depth temperatures are in excellent agreement with the 
measured values. 

The author wishes to express his appreciation to Barbara D. Arabian, 
Bette J. Stafford, and Davis D. Bland for their assistance in the pre- 
paration of the digital computer program. 


SYMBOLS 


0^ specific heat 

F exterior view factor 

F view factor-emissivity product to cabin environment 

env 

Hj heat of virgin, material degradation 

Hrp total enthalpy 

H wall enthalpy 

H^ q0 enthalpy of air at 3 00 ° K 

h film coefficient between backup materials 

h film coefficient between last backup material and cabin 

env ^ 

environment 

k thermal conductivity 


mass loss rate of char material 



i gas ablation rate 

g 

UP number of nodes in ablation material 

q hot wall convective heat flux with blowing 

C JDXOW 

q , heat flux due to combustion 

comb 

q cw cold wall convective heat flux without blowing 


^rad 

S 

S 

T 

T 

env 

r 

T 

T m 

VL 


AX 

A0 


e 

T1 

0 

e 


p 

a 


radiation heat flux 

surface recession depth 
surface recession rate 

temperature of node beginning of time step 
cabin environment temperature 

temperature of node at end of time step 
radiation heat sink temperature 
thickness of ablation material 
heat of combustion per unit weight of char 

thickness of a node 

time step (6' - 9) 

emissivity of material 

transpiration cooling efficiency 

initial time 

final time 

density 

Stephan-Boltzman constant 
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Subscripts: 

c charred state 

i node number 

j material number 

v virgin state 


PROGRAM DESCRIPTION 


The following general requirements were established before writing 
a digital computer program to analyze a charring ablation system: 

a. Stability of the equation for all applications. 

b. Machine running -time short enough to make use of the program 
economically feasible (a minimum of turn around time per problem). 

c. A minimum of input per problem. 

d. A wide variety of boundary conditions for application to both 
trajectory data and ground or flight test data analysis. 

STAB II has been formate ted in Fortran IV to analyze the transient 
thermal performance of a charring ablator heat protection system. The 
program considers one ablating material and up to 12 different backup 
materials with or without air gaps. Pure conduction or radiation and/or 
convection between backup materials is allowed. The ablation material 
may be divided into a maximum of 50 nodes and each backup material may 
be subdivided into a maximum of 10 nodes. The thermal properties of 
the materials are in tabular form and are temperature dependent. The 
ablation material is also dependent upon its state, that is, fully 
charred, partially charred, et cetera. 

The following surface boundary condition options are provided: 

a. Cold wall convective and radiative heat flux tables as a func- 
tion of time. These components are specified separately, since mass 
transfer at the surface blocks part of the convective heating, but in 
general, has no effect on the radiant heating. 

b. Surface temperature as a function of time. 
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c. Surface recession as a function of temperature or time. Slur- 
face recession as a function of temperature and pressure is also avail- 
able . 


Heat loss to the interior environment for the last node of the 
backup structure can be specified by two methods: 

a. Conduction into the node and radiation and/or convection loss 
to the interior environment. 

b. Conduction into the node and adiabatic wall. 


The STAB II numerical formulation of the equations describing the 
response of the heat shield is such that an implicit solution has been 
obtained. It is well known that numerical solutions of partial differ- 
ential equations are subject to several different types of errors. The 
first of these is the truncation error, due to the use of a finite sub- 
division. This error may be reduced by simply choosing a smaller sub- 
division, AX. The exact values are approached more and more closely 
as AX decreases. The second kind of error is the numerical, or round- 
off error. The way in which this numerical error grows or decays with 
time determines the stability of the difference equation. 


To illustrate the differences in the explicit and implicit equation 
form, consider a nonablating homogeneous solid. The one-dimensional 
Fourier conduction equation, neglecting any heat generation terms, is: 


SX 



( 1 ) 


The finite difference form of equation (l) written in the conventional 

,th 


forward time step or explicit form for the i node is 

( T i.i- T i) (H-Vi) r 

— “ = a5 


AX 


2k 


i-1 


AX 

2k i 


AX 

2k. 


AX 


P 


( 2 ) 


2k 


i+1 


where the prime superscript denotes values at .the end of the time step, 

A0 = e' - e 
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For explicit conduction solutions, the following stability criteria 
has been established: 


k 


A6 


^ 2 


which places an upper limit on the time step A 0 for a fixed truncation 
error. This criteria can require a prohibitive amount of machine time. 

Liebmann (ref. 3) advocated a solution of the equation which does 
not require this stability criteria. The finite difference equations 
are written in a backward time step form which affords an implicit so- 
lution. 

The implicit (backward time step) 
for the i node is: 

K-i - 0 « - 

AX’ AX_ “ AX 

2 k i-l 2k i 2k i 2 k i+l 

Equation (3) uses the temperature differences at the end of the finite 
time interval instead of the beginning, as in the explicit method, 
equation (2). T^ is the only known temperature in equation (3), but 

there are corresponding equations for each point in the system and all 
are solved simultaneously yielding the temperature at each node. 

Collecting all unknown temperatures on the left side of the equa- 
tion and the known temperature on the fight side, equation (3) is: 


difference form of equation (l) 


= Pc 


M K - T i) 

A6 


( 3 ) 



/ i + 1 

AX AX AX AX 
\ 2 k i-l 2k i 2k i 2 k i + l 


P 1 C P^\ , 

+ -a ir-J T i + 



P,c AX' 
1 P i 

“aF" 


(*) 
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Equation (4) is in the form of: 

AT^_ 1 + BT^ + CT^ +1 = D 


(5) 


STAB II generates such an equation for each node in the system. 

Since radiation is an important mode of heat transfer in charring 
ablative systems, a problem is encountered in any equation containing a 
radiation term. The radiation heat flux, written in a backward differ- 
ence form is: 

= Fea(r' ± k - tJ*) (6) 


This term cannot be used in an implicit solution since the unknown tem- 
perature T^ is to the 4th power. The 4th power unknown can be elim- 
inated by the following linearization: 


where 



(7) 


let 


v = ^ 


and rewrite equation ( 7 ) 


T i) 4 “ ( T i) U(1 + x)i 


If X has an absolute value near zero, the following is true: 

(1 + X) k as 1 + 4X 
Now substituting ( 9 ) into (8) 

( T i ) 4 * (h ) 4 ' 1 + - N'f + 4 f) ( T i f 

. 3 » 4 

s: 4T i ' J T i - 5T ± 


(8) 

(9) 


( 10 ) 
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Equation (10) is a linearized approximation of equation (7) in which 
the unknown temperature is only to the first power. The assumption in 
equation (10) is that has an absolute value near zero. Figure 1 

is a plot of the error obtained when (l + 4x) is substituted for 

(l + X) . For most ablation problems where the surface temperature is 
high and the radiation losses are significant, the value of T^ can 

easily be controlled to values of less than ±0. 1. 


Therefore, equation (6) can now be written 

^ “ Fea/W.V - JT. 4 - T^) 


(ID 


Using the linearized approximation for the radiation terms, the result- 
ing system of implicit difference equations constitute a tridiagonal 
matrix of the following form: 


B 1 T 1 + 


C 1 T 2 


D 


1 


A 2 T 1 + B 2 T 2 + C 2 T 3 


D 


2 


+ b ? t 5 + 


D, 


I 


A rf C lf-l + ®n t n " d n 

Gauss' elimination method, as discussed in reference 4, is applied to 
solve the system of equations. This method affords a fast and accurate 
solution for matrices containing a dominant diagonal. The solution of 
this matrix gives the temperature of each node in the system for the 
next future time step. The entire process is repeated for each time 
step throughout the run, giving a time history of the temperature at 
each node. 


Using this method of solution, residual errors in the temperature 
computations at the beginning of the time step are distributed through- 
out the entire system of nodal equations and tend to cancel out rapidly. 
The principle advantage in using the implicit method is a set of equa- 
tions that are mathematically stable in time and distance. Therefore, 
the magnitude of the time step is not limited by a convergence criteria. 
However, care must be taken in selecting the magnitude of the time step 
in order to minimize truncation errors when the second derivative of 
temperature with respect to time is large. A similar approach is used 
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to minimize truncation errors in distance toy choosing small node dimen- 
sions in locations where large second derivatives of temperature with 
respect to distance are expected. 

In the case of a char forming ablative heat shield where' approxi- 
mately 80 percent of the heat is reradiated, instability can arise in 
taking large time steps such that the temperature of the surface node 
can start oscillating on successive time steps on achieving a balance 
between the radiation source and sink. Therefore, in ablation problems 
in which the surface node loses a large percentage of heat by radiation, 
oscillations of the node can be damped out by taking small time steps 
during conditions of high heat flux and near radiation equilibrium tem- 
peratures. 


ANALYSIS 


Figure 2 is a schematic of the thermal protection system that is to 
be analyzed. A receding surface has been assumed with the formation of 
a residual char layer and reaction zone. The thermal protection system 
is composed of 1 charring material and a maximum of 12 different backup 
materials with or without air gaps. The analysis is such that the en- 
tire system may be composed of noncharring materials. The thermal prop- 
erties of all materials are temperature dependent; also, the charring 
material properties are state dependent (fully or partially charred). 

The response of charring ablation heat shields to a hyperthermal 
environment is extremely complex, and simplifying assumptions and approx- 
imations are necessary to afford a numerical solution. The following 
assumptions and approx ima tions are utilized in the equations developed 
in this report: 

a. The material decomposes from the virgin state to a porous char 
layer in the reaction zone. 

b. The reaction zone can be defined by an upper and lower temper- 
ature limit. 

c. The gas generated within the reaction zone is assumed to pass 
out of the structure with no pressure loss. No gas accumulation within 
a node is allowed. 

d. Local thermal equilibrium is maintained between the gas and 
solid. 

e. The gas undergoes no further chemical reaction within the re- 
sidual material after having been formed. 
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Derivation of Equations 


The equations are derived for a moving "boundary coordinate system, 
where the front face is the moving surface (ref. 5). With this system, 
the ablating material is divided into a fixed number of nodes with a 
thickness (AX) which depends on the instantaneous location of the 
front face. The surface recession is handled in a continuous manner 
eliminating the need of throwing away or lumping off of nodes. 

The physical model for the front surface including all heating 
terns is shown below 


AX 

2.0 



The energy equation at the front char surface is: 

dT, 


A ' “ Pl Vl) • * * ^ + ‘ ^ 


= Q. + m c T 0 + p 0 S_c T_ - m c T. 

m g 2 p 2 2 r 2 2 p 2 2 g x p x 1 


‘ P l C Pl \ T l " ^l-2^Ax) 


( 12 ) 


where 


^in _ ^c Blow + < ^rad + ^comb ” Fecr (^1 “ ^°° ) 
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and 


a(AX) a ( yl - s \ _ s 

ae a5 \rp - ly “ rp - 1 


•where S is the linear surface recession rate, NP the total number 
of nodes in the ablation material of thickness VL. 

Rewriting equation (12) in implicit finite difference form: 

. K - T j) 

Q. +m c T_ - m c T. - Sp.c T. - 7= rr- 

in gp pp 2 g- P- 1 1 P- 1 

2k 1 2k 2 


p e g m - ^ T - 

P 2 p 2 ^RP - 1.0 J X 2 


P-, c 


m - T 
AXl 1 X 1 


lp 1 2 \ W 


Rearranging and collecting terms: 


, • AX 

-I me + Sp., c + P_c T75 + 7=7 

I g-L P-l 1 P x 1 p x 2A6 AX 


- (ipTl) ^ 

1 .ipe f— M\ t' 

+ £X_ 2 p lp 1 \RP - V I 


2k^ 2kg 


+ ' V% + STT^I + p 2 c p/ (® 

/*Y»_ * 


RP - 


1 il\\ T ’ 


l.oy 2 


2k 1 2k 2 


£X 


= - p lV 2A5 T 1 - Q 


in 


(12b) 
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The physical model for interior points in the mature char zone, in- 
cluding all heating terms, is shown below: 



The energy equation for interior points in the char matrix is: 


d dT. t / • \ 

35 ’ Mp i c Pl a® - - h c Pi T i (bpTt) 


“A c T. _ + k 

g i+l p i + l 1+1 


( £P\ 

i-1, i VAX/ 


UP - i - | , 

+ P i+l C p i+1 S V UP - 1 / T i+1 


, / AT\ . ' t ( KP - i + ^ \ , 

k i,i+l \AX/ ' m g. C p. T i - P iV. S \ NP - 1 ) T i 


3 i *1 


( 13 ) 


IrolH 
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Putting equation ( 13 ) in an implicit finite difference form: 


/ 1 

AX 

u 42L 

\ 2k i-l 

2k i 


. / . [HP - i + % 

T, - /&. e + P^S 



AX + AX 


2k i _ 1 2k ± 


1 AX 

ax . 5 ~ + p iV 35 • p iV *1 

2k i a i + i 


BP - i - 


+ P-.^c 


iAl %- + #- i+1 p i+ i V ^ " 1 


2k i 2k i + l 



AX 


= “ P i c p. Z5 


T 

P i AO i 


(13a) 

NOTE: In the mature char zone, no internal gaseous ablation products 

are assumed to form. The reaction zone is the source for the formation 
of the. internal gaseous products. Therefore, in equations (12). and (13), 
m = m 
g i g i+l 

The physical model for nodes in the reaction zone is identical to 
schematic shown for the interior nodes in the char except for consider- 
ing the energy absorbed in formation of the gaseous ablation products. 

The heat balance equation for a node in the reaction zone is: 


cLT / • \ 

s^V 1 ) = ^ ' Pi Vi ' {\ ' \ +1 ) 




= m c T. ,. + k 
®i+i *1+1 i+1 


/ AT\ 
1-1,1 y AXy 


NP - i - f , 

+ P i+l C p i+1 S V HP - 1 / T i+1 


_ k (££ \ m i c T* - P.c S 

k i,i+l \AX/ % Pi 1 1 p i 


r NP - i + - , , 


V 

2 


BP - 1 


(iM 



Ik 


rearranging: 




P,c S 

1 p i 



+ 


1 


AX 

2k i-l 



2k i 


AX 


2k 


n & 

p i c Pi a5 - 


i+1 



T i + 



+ 


1 



AX 


2k 


i+1 


+ 





- m 
g 


H a 


i+1, 


(l!+a) 


The physical model for the interface between the reaction zone and vir- 
gin material is illustrated "below: 


— AX — ► 
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The heat balance equation for this node is: 


dT / • \ 

OT (®iVi) ■ ^V Pl it - V Pl *i (iipTT) - 


\ /BP - i - f ] , 

- (I) + Vi°p 1+1 s U-iV 


- k (^\- 6 c t' 
k i,l+l\£KJ g« p. i 


J w - i + t |„- 

P i C Pi S \ BP - 1 J T i 


Rearranging: 


(15) 


AX AX 1 i-1 


T. , - m c + 


1 + — i 


,^1 2k. 


g, p, + 

^i-i 2k i a i 2k i+i 


UP - i + 2 
+ P i C p J* \ BP - 1 


J A AX 

+ P iV a5 


- p i c p, (bp - lj) 1 


1 

AX AX + P i+l C p S \ BP - 1 
■ 2k i 2k i + l 


* T i + l 


AX m . TT 
= -P.c rs T. - m H , 
l Pi Ao l g i d 


(15a) 
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The physical model for an interior node in the virgin material is; 

AX— — i 


S. . .c P.T."* 
i-l,x p. i i 


/at\ 
C i-1, i \AXy 



The heat balance for this nonablating node is: 


, . dT ' , . A \ 

d5 («V Vi) - dT - ’i'p/i (spvr) 


= k 


/AT N 

L-l,i VAX/ + P i+l C p^., S \ HP 


. f® - i-i\ . 


i+1 


112 • (M\ 

-1 / i+1 *i,i+l\AX/ 


'HP - 


- p i C Pi S \— 




\ (16) 



IT 


Rearranging: 



(l6a) 


The physical model for the last node in the ablation material and first 
node in the backup structure is: 
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The heat "balance equation for this node: 

AX 


aMh* + °Pi, 3+1 P i.J + ii T iJ 


/'AX.p. .c + AX.,.c p. ..A 

J i,J P i<t1 j +1 Pj.j-n i ^ +1 ) dT 1 / S_\ » 

V 2 / dS “ 2 Up - i / c P 1# j p i,3 T i 


= k 


AT 

1-1,3 

1,3 


1 

' A, 3 S \w^rjJ T i 


- k 


1,3+1 

i+1, 3+1 


/at'. 


(17) 


Rearranging: 


T. 


AX. AX. \ "i-l 

^ + 


2k . _ . 2k . . 
i-i* J i,3 


AX. AX. 

, il 4. ■ J 

2k. , . 2kj 


1-1,3 1,3 


fa 

2k i,3+i ^1+1,3+! 


/AX.P. c + AX c f 

2 J+1 . p j ,l+i 


i,3+l\ i , 

T i + 




)‘ I+1 


2k i,3+l 2k i+l, 3+1 / 


= _ ^ j A ./ 1 * 3 J+lCp 1 .. 1 +l Pl,3+1 j T> 


(17a) 

The backup structure may contain up to a maximum of 12 different 
materials with or without air gaps between materials. Therefore, con- 
duction or radiation and/or convection between materials is allowed. 

The heat balance equations for the various modes of heat transfer in the back 
backup structure are now presented: 
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a. Interior node material: 



(18a) 

t. First node of two interior materials with no gap: 



NOTE: 


T = T 


for this case 
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Rearranging: 



Equation (20) may "be linearized using the approximation 

T* a: 4 tV - 

as discussed in the Program Description section. 
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Therefore, rearranging and linearizing, equation (20) "bee ernes: 




d. Last node of an interior material with an air gap between 
materials: 






( 21 ) 
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Rearranging and linearizing: 




e. Final node in backup structure: 
(l) Adiabatic surface 



Rearranging: 




(22a) 
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(2) Radiation and/or convection loss to cabin environment 



Rearranging: 



(23a) 

Discussion of Assumptions 

A brief discussion of several assumptions and approximations made 
in deriving the heat balance equations is now presented. 

As shown in the Derivation of Equations section, transient heat 
conduction, thermal degradation, and the flow of the gaseous products 
from the reaction zone are the internal thermal transport phenomena of 
interest. Several methods are available in the treatment of the thermal 
decomposition process and they differ primarily in whether the chemical 
decomposition occurs in a single plane at a fixed temperature or a 
spacially continuous decomposition in depth is assumed. This analysis 
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assumes that the decomposition from the virgin to the char state occurs 
in a reaction zone that is defined by known temperature limits. These 
temperature limits are determined from the rmogravime trie (TGA) test 
data for the particular material being investigated. Figure 3 is a TGA 
curve for typical charring ablation material. From this curve, the 
rate of pyrolysis is calculated by knowing the temperature change 

of a particular node with time, that is. 


p i = 


P i " P i 

“SB — 


(24) 


HP 

m = 2 p AX (25) 

s i i 1 


This method of computing the gas generation rates and local instantan- 
eous density may be subject to error since the TGA curve of a material 
is influenced by temperature rise rate (DEG/SEC) and the reaction zone 
may shift up and down the temperature scale. This error can be elim- 
inated by the use of an Arrhenius expression of the foxra 

|| . -A(p - p c ) n p HT (26) 

The method now being used in STAB II (equation (25)) to calculate the 
pyrolysis rate is being investigated to determine its validity. The 
final formulation of the pyrolysis rate law must rest heavily on the 
experimental rate data for the material under investigation. The use 
of simple expressions such as equations (24) and (25) may be entirely 
adequate depending upon activation energy for the decomposition pro- 
cess and order of reaction. 

The aerodynamic heating input in the analysis consists of convec- 
tive and radiative components treated separately. This distinction is 
necessary since the convective heating can be significantly reduced 
due to the injection of the ablation gases into the boundary layer with 
generally no effect on radiant heating. Reduction in the convective 
heating rate can be approximated by the following expression (ref. 6) 

W s Wgfa - H v) 


(27) 
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Therefore, 


Sp - H 


%lo» - 4 ow ^ 


- q. 


Block 


(28) 


However, equation (28) is unsatisfactory for high blowing rates, since 
%lock can ^ ecame greater than 4 CW * An experimental curve of block- 
ing effectiveness J as a function of the mass transfer param- 

\*c vj 

eter v — can be employed to deteimine the heating reduction at high 
hew 

blowing rates. Both methods have been employed in the STAB II analysis. 
Equation (28) is presently in use. However, no satisfactory method for 
accurately predicting the convective heat blockage has been determined. 

Another source of heating is the combustion of the ablation prod- 
ucts in the boundary layer. Reference 7 presents an analysis of the 
oxidation of a carbon surface and the resulting combustive heating. The 
heating due to combustion as derived in reference 7 is 


^■camb ~ “c ^c 


(29) 


where £H c is the heat of combustion per unit weight of char. 

The thermal properties of the ablation material are both tempera- 
ture and state dependent (fully or partially charred). Figure 4 is an 
illustration of the variation of these properties with temperature and 
state. The thermal properties are assumed to vary as follows: • 

a. Char zone (r ± > 

k Q = f (temp) 

c = f (temp) 
p c 

p = constant 
c 

b. Reaction zone ^T^ * T. < T^j 

p = f (temp) = P v + (P v - P c ) 
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c. Virgin zone j 

p = constant 
v 

k y - f (temp) 

c = f (temp) 
p v 

The calculation of char remo-val, due to chemical, thermal, mechanical, 
or "by comhination of these mechanisms, has "been examined by a multitude 
of investigators and numerous correlations exist, depending on the 
specific material involved. 

To provide a maximum degree of flexibility for analyzing both 
ground and flight test data and systhesizing trajectories, the follow- 
ing provisions for char removal (surface movement) are provided: 

a. Removal of char as a function of surface temperature. 

b. Removal of char at a rate which is a function of time. 

As the char is removed, the surface moves with respect to a coordinate 
fixed in the material. The distance between the initial surface loca- 
tion and the char surface is 

r e • 
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CUSTOMER UTILIZATION INSTRUCTIONS 


Introduction 

IBM 7094/40 program 1021, standard ablation program, designated 
STAB II, is designed to evaluate the transient thermal performance of 
a charring ablation heat protection system. The program considers 1 
ablating material and up to 12 different materials in the supporting 
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backup structure. A maximum of 50 nodes may be considered in the abla- 
tion material and a maximum of 10 nodes per material is allowed for 
each backup structure material. Air gaps can be considered between 
successive materials in the backup, thus allowing for both radiative 
and/or convective heat transfer between materials. The heat loss to 
the cabin environment from the backup structure can be accomplished by 
both radiation and/or convection or an adiabatic backface surface may 
be prescribed. 

Unless otherwise specified, the input problem data is in "floating 
point" form (E12.8 format) and must end in columns 12, 24, 36, 48, 60, 
and 72. It is suggested that each floating point number have a sign, 
a two-digit exponent, and a decimal point. For example, the number 
145.23 can be written as +1.4523 +02, +145.23 +00, or +14523 +03 

Input Nomenclature 

The nomenclature used in the problem data input is as follows: 

NCASE number of problems to be run successively 

HEAD any 72 alphabetical and/or numerical characters 

TITLE control card for reading in new input for successive 
problems 

1. blank card — new data will be read in 

2. 6 asterisks in columns 1 to 6. Skip to next 

read statement 

TLIM time limit of problem, sec 

TINT starting time of problem, sec 

NPTT number of points in time-step table (the minimum 

value of NPTT is 2) 

NPL0T output plot control 

=1 plot routine will be used 

=0 plot routine will be ignored 

TTABLE time in time-step table, sec 
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DELTT 

IPRC 

PC0EV 

ERAD 

TABL 

TCHAR 

TREC 

RHj0V 

RH0C 

EBL0W 

EMV 

EMC 

H500 

YL 

HY 

YPT 

=0 

=1 

W 


time step to be used for each calculation — starting 
at time TTABEE, sec 

variable print frequency in TTABEE table; that is, if 
DELTT = 1. 0 and IPRC =10, the output will be printed 
at 10- second intervals 

factor to correct convective heating rate for various 
body locations 

factor to correct radiative heating for various body 
locations 

temperature at which ablation starts, °R 

temperature at which ablation stops, °R 

surface temperature, (°R) or time (second) at which char 
removal is to start 

density of virgin ablation material, lb / ft'’ 

density of mature char material, lb/ft^ 

blowing efficiency of ablation gases in reducing convec- 
tive heating 

emissivity of virgin ablation material 

emissivity of charred ablation material 

enthalpy of air at 300° K, 129*06 Btu/lb m 

initial thickness of virgin ablation material, in. 

heat of degradation of virgin material, Btu/lb 

test to determine if the reaction zone and char zone 
thermal properties are irreversible with temperature 

properties are irreversible and equal to the value at 
the maximum individual node temperature (this is the 
recommended value for VPT) 

properties are reversible 

view factor for external environment 
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TV 

CHARK 

CHARC 

ABLK 

ABLC 

HP 

NKC 

NCPC 

HECV 

NCPV 

HREC 

TKC 

XKC 

TCPC 

CPC 

TKV 

XECV 


sink temperature of external environment, °R 

thermal conductivity of material at TCHAR, Btu/hr-ft-°R 

specific heat of material at TCHAR, Btu/lb m -°R 

thermal conductivity of material at TABL, Btu/hr-ft-°R 

specific heat of material at TABL, Btu/lb m -°R 

number of node points in ablation material 

number of points in char thermal conductivity — temper- 
ature table 

number of points in char specific heat — temperature table 

number of points in virgin thermal conductivity — 
temperature table 

number of points in virgin specific heat - temperature 
table 

number of points in surface recession — temperature or 
time table 

temperature values in char thermal conductivity — 
temperature table , °R 

thermal conductivity values in char thermal conductivity — 
temperature ' table , Btu/ ft-hr- °R 

temperature values in char specific heat - temperature 
table, °R 

specific heat values in char specific heat — temperature 
table, Btu/lb m ~°R 

temperature values in virgin thermal conductivity — 
temperature table, °R 

thermal conductivity values in virgin thermal conduc- 
tivity — temperature table, Btu/ft-hr-°R 

temperature values in virgin specific heat — temperature 
table, °R 


TCFV 
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CPV 

TS 

SR 

ETRAPT 

TIME 

QRAD 

VEL 

KMB 

EPBS 

BL 

XKPM 

NKPB 

NCPB 

XIDKT 

TXK 

XK 


specific heat values in virgin specific heat temperature 
table, Btu/lh m -°R 

temperature (°R) or time (sec) values in the surface re- 
cession table 

surface recession values in the surface recession — 
temperature or time table, in. /sec 

number of time points in the trajectory input table 

the array of (ETRAPT) trajectory time values, sec 

the corresponding array of cold wall convective heating 
2 

rates, Btu/ft -sec 

the corresponding array of radiative heating rates, 

O 

Btu/ft -sec 

the corresponding array of flight velocity, ft /sec 

number of materials in backup structure 

total number of node points in backup structure 

total thickness of backup structure, in. 

number of nodes in each individual material in backup 
structure 

number of points in each individual backup structure 
material thermal conductivity — temperature table 

number of points in each individual backup structure 
material specific heat — temperature table 

any 72 alphanumeric characters used to describe each 
individual material in the backup structure 

temperature values in backup material thermal conduc- 
tivity — temperature tables, °R 

thermal conductivity values in backup material thermal 
conductivity — temperature tables, Btu/ft-hr-°R 
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TCP 

CPX 

RH0BX 

XBM 

EMFB 

EMBB 

H 

GAPX 

PTEST, 

BTEST 

=0 
=+l 
— 1 
TENV 
HENV 

EEKV 

QL0SS 

=0 

=+l 


temperature values in Backup material specific heat - 
temperature tables, °R 

specific heat values in Backup material specific heat - 
temperature tables, Btu/lB^- °R 

"Z 

density of individual materials in Backup, lb/ftr 

thickness of individual materials in Backup, in. 

emissivity of front surface of each material in Backup 

emissivity of Back surface of each material in Backup 

film coefficient Between adjacent materials in Backup, 
Btu/hr-ft 2 -°R 

width of gap Between adjacent materials in Backup, in. 

tests to determine the mode of heat transfer Between 
materials for the front and Backface of each .material 
respectively 

conduction only Between materials 
convective heat transfer only 

radiation only or radiation and convection heat transfer 

temperature of interior cabin environment, °R 

film coefficient to interior cabin environment, 
Btu/ft 2 -hr-°R 

view factor and emissivity product for radiative heat 
transfer to cabin interior 

Boundary condition Between last node of the Backup 
structure and cabin environment 

adiabatic surfaces 

radiation and/or convective loss 
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TEST2 determines the proper heat shield initial temperature 
distribution 

=0 constant, uniform initial temperature distribution 

=-l arbitrary initial temperature distribution 

=+l linear temperature distribution 

TEMPI temperature to be used when constant temperature dis- 
tribution option is used, °R 

TX0 initial temperature at front surface of heat shield to 

be used in computing initial linear temperature 
gradient, °R 

TEMPI arbitrary temperature distribution values, to be used 
only if TEST2 is negative, R 

NHP number of points in enthalpy — temperature curve fit 

HX enthalpy values in enthalpy — temperature table, 

Btu/lb 

TW corresponding temperature values in enthalpy — tempera- 

ture table, °R 


NOTE: This table is used for computing the wall enthalpy. 

An input deck for NHP, HX, and TW has been prepared for air and is 
available upon request. 


Input Data Card Preparation 

The input data are given in the following order. Each number 
below refers to a separate record and must begin on a new data card. 

The input data has been grouped, where possible, into various sections 
dealing with a particular part of the input, that is, ablation material 
properties, trajectory data, backup structure, et cetera. This permits 
the use of a minimum number of input cards for running successive prob- 
lems. The title card as described in the input nomenclature controls 
the input for successive problems. 

1. The first data card contains the value of NCASE. NCASE is 
an integer (15 format) and must end in column 5* This card tells how 
many problems are to be run and is entered only once at the start of 
the data deck. 
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2. Columns 1 through 72 of the second data card contain any title 
or identification information desired — any alphanumeric character may 
he used. This card vill he printed at the top of the first page of the 
output. This card must he included in all successive problems to he 
run. 


a. Problem Time Section 

5. TITLE card — if blank, the following 2 cards must he submitted; 
if 6 asterisks are punched in columns 1 through 6, skip to record num- 
ber 6. 

4. This record contains, in the following order: TLIM, TIM 1 , 

NPTT, NPL0T. TLIM and TIM are entered as floating point numbers and 
must end in columns 12 and 24. NPTT and NPL0T are integers entered 
with an 15 format and must end in columns 5° and 55* 

5. Start entering the values of TTABLE, DELTT, IFRC. TTABLE and 
DELTT are floating point numbers and must end in columns 12 and 24. 

IPRC is entered as integer with an 15 format and must end in column 5 0. 
Use as many cards as required to enter NPTT values. 


b. Heating Rate Factors Section 

6. TITLE card — if blank, the following card must be submitted; 
if 6 asterisks are punched in columns 1 through 6, skip to record num- 
ber 8. 

7. Enter the following: PC fiN, FRAU. These numbers are entered 
as floating point numbers and must end in columns 12 and 2k. 


c. Ablation Material Section 

8. TITLE card — if blank, the following cards 9 through 18 must 
be submitted; if 6 asterisks are punched in columns 1 through 6, skip 
to record number 19 * 

9- HEADNG card — any alphanumeric characters in columns 1 through 
72. Records 9 through 18 contain input data for the ablation material. 

10. Enter the following: TABL, TCHAR, TREC, RH0V, RH0C, FBL0W. 
These numbers are entered as floating point numbers (6E12.8 format) 
must end in columns 12, 2k, 5 6 , 48, 60 and 72. 
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11. Enter the following: EMV, EMC, H3OO, VL, HV, VPT. Use same 

format as card 10. 

12. Enter the following: FV, TV, CHARK, CHARC, ABLK, ABLC. Use 

same format as card 10. 

13. This card contains, in the following order: UP, NKC, NCPC, 

NKV, NCPV, NREC. These numbers are fixed point integers and must end 
in columns 5, 10, 15, 20, 25, and 30. An 15 format is used to read in 
these numbers. 

14. Start entering the curve of TKC versus XKC, with the values of 
TKC ending in columns 12, 5 6, and 60. The corresponding values of XKC 
must end in columns 24, 48, and 72; for example, three TKC -XKC points are 
contained on one card. The numbers are entered as floating point num- 
bers. Use as many cards as required to enter NKC points on the curve. 

15. Start entering the curve of TCPC versus CPC with the values of 

TCPC, ending in columns 12, 36, and 60. The corresponding values of 

CPC must end in columns 24, 48, and 72; for example, three TCPC-CPC 

points are contained on one card. The numbers are entered as floating 
point numbers. Use as many cards as required to enter NCPC points on 
the curve. 

16. Start entering the curve of TKV versus XKV with the values of 
TKV ending in columns 12, 36, and 60. The corresponding values of XKV 
must end in columns 24, 48, and 72; for example, three TKV-XKV points 
are contained on one card. The numbers are entered as floating point. 

Use as many cards as required to enter the NKV points on the curve. 

17. Start entering the curve of TCFV versus CFV with the values of 

TCPV, ending in columns 12, 36, and 60. The corresponding values of 

CFV must end in columns 24, 48, and 72; for example, three TCFV-CFV 

■points are contained on one card. The numbers are entered as floating 
point. Use as many cards as required to enter NCFV points on the curve. 

18. Start entering the curve of TS versus SR with the values of TS, 
ending in columns 12, J>6, and 60. The corresponding values of SR must 
end in columns 24, 48, and 72; for example, three TS-SR points are con- 
tained on one card. The numbers are entered as floating point. Use as 
many cards as required to enter NREC points on the curve. 


d. Trajectory Data Section 

19. TITLE card — if blank, the following cards 20 through 22 must 
be submitted; if 6 asterisks are punched in columns 1 through 6, skip 
to record number 23. 
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20. HEADKG card — any alphanumeric characters in columns 1 through 
72. Records 21 through 22 contain trajectory input data. 

21. Enter the following: KTRAPT. This number is an integer and 
must end in column 5- An 15 format is used to read in this number. 

22. Start entering the trajectory data in the following order: 
TIME, QC0N, QRAD, VEL. These values are entered as floating point num- 
bers and must end in columns 12, 24, J>6, and 48. There are four tra- 
jectory data points on one card. Use as many cards as required to 
enter NTRAPT points in the trajectory. 


e. Backup Structure Section 

25. TITLE card — if blank, the following cards 2k through J>2 must 
be submitted; if 6 asterisks are punched in columns 1 through 6, skip 
to record number 53* 

24. HEADKG card — any alphanumeric characters in columns 1 through 
72. Records 25 through 52 contain properties of backup structure. 

25. Enter the following: KMB, HPBS, BL. These three values must 

end in columns 5* 10 and 24. KMB and KPBS are integers and are read in 
under an 15 format. BL is a floating point number. 

26. Enter the values of XKPM. XKBi is in floating point form and 
must end in columns 12, 24, 56, 48, 60 and 72. Use as many cards as 
required to enter KMB points. 

27. Enter the values of KKPB and KCPB. These numbers are integers 

and KKPB must end in columns 5* 15 > 25, 55* and 45* and the correspond- 
ing values of KCPB must end in columns 10, 20, 40, and 50. An 15 

format is used to read these values. Five KKPB-KCPB values are con- 
tained on one card. Use as many cards as are required to enter KMB 
points. 

28. XIDHT card — any alphanumeric characters in columns 1 through 
72. This card contains a description of each backup material. 

29. Start entering the curve of TXK versus XK with the values of 
TXK, ending in columns 12, 56, and 60. The corresponding values of XK 
must end in columns 24, 48, and 72; for example, three TXK-XK points 
are contained on one card. The numbers are entered as floating point. 
Use as many cards as required to enter KKPB points on the curve. 
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30 . Start entering the curve of TCP versus CPX with the values of 
TCP, ending in columns 12, 36 , and 60. The corresponding values of CPX 
must end in columns 24, 48, and 72; for example, three TCP-CPX points 
are contained on one card. The numbers are entered as floating point. 
Use as many cards as required to enter NCPB points on the curve. 

NOTE: Repeat records 28 , 29, and 30 until the properties for MB 

materials have been entered. The maximum number for MB is 12. 

31. Start entering the following values in order: RH0BX, XBM, 

EMFB, EMBB. These values are entered as floating point numbers 
(6E12.8 format) and must end in columns 12, 24, 36 , 48, 60, and- 72. Use 
as many cards as required to enter MB points. 

32. Start entering the following values in order: H, GAPX, FIEST, 

BTEST. These values are entered as floating point numbers (6E12.8 for- 
mat) and must end in columns 12, 24, 36 , 48, 60, and 72. Use as many 
cards as required to enter MB points. 


f. Interior Environment Section 

33 . TITLE card — if blank, the following cards 34 through 35 must 
be submitted; if 6 asterisks are punched in columns 1 through 6 , skip 
to record number 36 . 

34. HEADNG card — any alphanumeric characters in columns 1 through 
72. Record 35 contains properties of environment. 

35* Enter the following: TEW, HENV, FEW, QL0SS. The values are 

entered as floating point numbers and must end in columns 12 , 24, 36 , 
and 48. 


g. Initial Temperature Section 

36 . TITLE card - if blank, the following records 37 through 38 must 
be submitted; if 6 asterisks are punched in columns 1 through 6 , skip to 
record 40. 

37. HEADNG card — any alphanumeric" characters in columns 1 through 
72. Records 38 through 39 contain initial temperature distribution in- 
put. 


38 . Enter the following: TEST2, TEMPI, TXO. These values are 

entered as floating point numbers and must end in columns 12 , 24, and 

36. 
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NOTE: If TEST2 is a negative number, record 39 must be submitted* 

otherwise, skip to record 40. 

39- Enter the arbitrary temperature distribution values, TEMDI. 
These values are entered as floating points with a 6E12.8 format. Use 
as many cards as required to enter NP plus NPBS node points. 


h. Enthalpy — Temperature Section 

40. TITLE card — if blank, the following records 4l and 42 must be 
submitted* if 6 asterisks are punched in columns 1 through 6, this is 
the last data card in the problem input. 

41. Enter the following: NHP. This value is an integer and must 

end in column 5* An 15 format is used to read in this number. 

42. Start entering the curve of HX versus TW with the value of HX 
ending in columns 12, J>6, and 60. The corresponding values of TW must 
end in columns 24, 48, and 72; for example, three HX-TW points are con- 
tained on one card. The numbers are entered as floating point. Use as 
many cards as required to enter NHP points on the curve. Record 42 
consists of the last data cards required as input for a problem. 

As many successive problems as you wish may be run at one time by 
proper input preparation. STAB II has been designed to save all input 
information until changed by new input data. Therefore, the use of the 
TITLE control card is very important .when running more than one problem 
and using the input data of the previous problem(s). As shown, each 
input section starts with a TITLE control card for determining whether 
new input data is to be used. ■ If any data is changed within a section, 
then all data cards required for that section must be submitted. 

■ STAB II can also be used for solving one dimensional transient heat 
conduction problems of nonablating materials. The following input param- 
eters must be adhered to: 

a. TABL must be greater than the maximum temperature expected dur- 
ing the calculation. Also, TABL > TCHAR > TREC. 

b. The ablation material must be considered to be the first mater- 
ial in the structure for calculational purposes. 

c. The virgin and char properties must be inputed as described 
above but can have the same values; that is, XKV = XKC, CPC = CFV, 

RH0V = RH0C, et cetera. 
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The following dimensions statements and program limitations should 
not he violated when preparing the input described above for ablating 
and nonablating structure: 

a. All property tables can have a maximum of 20 points (i.e., a 
temperature and specific heat value constitute one point). 

b. The surface recession table can have a maximum of 50 points 
(TS and SR constitute one point). 

c. The trajectory table can have a maximum of 3 00 points (TIME, 
QC0N, QRAD, VEL constitute one point). 

d The ablation material can be broken into a maximum of 50 nodes. 
The backup structure can consist of up to 12 different materials with a 
maximum of 10 nodes per material. 

e. A minimum of 3 nodes per material (ablation or backup) must be 
specified.. 

f. A minimum of two materials must be specified (ablation material 
and one backup structure material). 

g. Pure conduction only is allowed between the ablation material 
and the first material in the backup. 

h. If any data input is changed in the Ablation Material Section 
on successive problems, the Ablation .Material Section data cards plus 
the Initial Temperature Section data cards must be submitted. 


Program Output Information 


The computed results are available in two forms of output] tabular 
and plot outputs. The tabular output presents the computed results in 
block type form for each computation step as controlled by the print 
count control number. As discussed in the preparation of input data, 
both the computational time step and print control can be varied through- 
out the running of a problem. Therefore, excessive printed output is 
avoided, as well as a considerable savings in actual machine computa- 
tion time. The plot outputs are printed and plotted only when the en- 
tire set of problems to be run are completed. 
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Tabular output . - The program prints a listing of the data input 
parameters for identification of the problem and ease in determining if 
there are any input mistakes. For stacked problems, the program prints 
only that input information that is changed from the previous problem. 
The following calculated problem output is printed: 

a. Time, sec 

2 

b. Cold wall convective heating rate without blowing, Btu/ft -sec 

2 

c. Radiative heating rate, Btu/ft -sec 

d. Velocity, ft/sec 

2 

e. Gas ablation rate, lb /ft -hr 

HI' 

f. Char ablation rate, lb /ft 2 -hr 

Hr 

2 

g. Total ablation rate, lb. /ft -hr 

Hr 

h. Surface recession depth from original surface, in. 

i. Hot wall convective heating rate without blowing, Btu/ft -sec 

j. Temperature distribution in ablation material, °R 

k. Temperature distribution in backup structure, °R 

The temperatures printed for the ablation material are for fixed 
distances from the original surface. These distances are calculated 
from the initial ablation material thickness and number of nodes in 
ablation material. For example: 

let 


VL = 1. 0 in. 


HP - 11 


then 


AX = 


VL 

HP - 1 


0.1 


The temperatures will be printed for X distances of 0, 0.1, 0.2, 
0.5, et cetera, from the original surface until the surface has receded 
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"beyond these fixed distances at which time the node no longer exists 
and is dropped from the printout. This is illustrated in the following 
way using the above example} let surface recession = 0.26 in then the 
first temperature printed is the surface temperature of the material, 
located 0.26 in. from the original material surface. The following 
printed ablation material temperatures are for X distances of 0.3, 0.4, 

0 . 5 , .... 1.0 in. 

The for ma t for the temperature distribution printout is E16.5 with 
six temperatures printed per line. 

Plot output . - The plot output gives the following ablative material 
performance parameters as a function of time: 

a. Surface depth, in. 

b. Bondline temperature between ablator and backup structure, °R 

c. Two selected isotherm depths 

These values are also printed in tabular form for ease in checking 
and replotting of the results. The plotted curves contain all maximum 
and minimum values of the parameters. 


PROGRAM •VERIFICATION 


As discussed in the previous sections, approximations and assump- 
tions were made in the analytical model to afford both a quick and 
accurate solution in predicting the thermal response of a charring heat 
shield. These simplifying assumptions and approximations are expected 
to introduce only minor errors; however, the validity of the analyses 
and resultant accuracy can be judged only by a comparison with exact 
theoretical solutions and experimental data. Three examples have been 
selected and a comparison of the STAB II results with the theoretical 
and test data is discussed in the following paragraphs. 

An elementary transient heating example was chosen to demonstrate 
the accuracy and numerical stability of the STAB II program. A steel 
slab 6 inches thick was selected and assumed to be at uniform initial 
temperature of 460° R (0° F). The thermal properties were considered 
constant. The front surface was subjected to a heating rate of 

72 Btu/sec-ft and an adiabatic back surface was assumed. Figure 5 
shows a comparison of the STAB II calculated in depth temperatures as 
a function of time with the exact solution taken frcm reference 8. 
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To demonstrate the STAB II solution with a moving boundary, a slab 
with constant properties, uniform initial temperature, front surface 
moving with a constant velocity, and constant surface temperature 
was chosen. The exact solution for semi-infinite slab with these bound- 
ary and initial conditions is presented in reference 9* ' Figure 6 pre- 
sents a comparison of the STAB II temperature response with the exact 
solutions. As can be seen from this figure, the two solutions are not 
in agreement for approximately the first 50 to 60 seconds of the trans- 
ient. This disagreement is the result of the quasi-steady state 

((l) 5=0 ■ °> k (S) x=0 ’ ® PC P ") 

assumption made in the exact solution analysis. A calculation was made 
to estimate the induction time (time at which ^ = 0 is a good assump- 
tion) and found to be approximately 60 seconds, which is in agreement 
with the STAB II results. 

Finally, to verify the fully charring ablation model, an example 
of a typical charring material was chosen. The charring ablation mater- 
ial is initially 1.6 inches thick with an adiabatic back surface and a 

2 

constant heat flux of 95 Btu/sec-ft applied to the front surface. The 

surface is assumed to recede at a constant velocity of 3*05 (l0 - ^) in. /sec. 
Figure 7 presents a comparison of the in-depth temperatures with actual 
test results obtained in an arc tunnel. The results are in good agree- 
ment, with the largest deviations between calculated and measured values 
occurring for the thermocouple located 1.0 inches in depth. The dis- 
agreement could be attributed to several possible errors; thermal prop- 
erty values, incorrect location of thermocouples, et cetera. The effect 
of varying the thermal properties (thermal conductivity, specific heat, 
etc. ) is presently being investigated and will be reported in a future 
report. 

Tables I and II present the input and output data used for this 
example. Figures 8, 9> and 10 are the resulting plot routine output. 

The camparisons presented above between the computer results and 
the exact solutions and test results are considered satis factory. As 
discussed previously, the assumptions used in the analytical model will 
be examined more critically as additional test data and analyses become 
available. The program will be revised and updated as required to re- 
flect these additional studies. 
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An analyses and computer program for predicting the transient 
thermal response of a charring ablation thermal protection system has 
been described. The numerical formulation of the equations is such 
that an implicit solution is obtained. This method of solution affords 
both a rapid and accurate solution for both ablating and nonablating 
type problems. 

Provision is made in the program for a number of surface boundary 
conditions. These provisions allow efficient use of the program for 
analyzing both ground and flight test data and trajectory synthesis. 

The computer program has been checked out with both exact solu- 
tions and actual ablation test data. The numerical results are in good 
agreement with the exact solutions and test data. However, the analysis 
depends upon using good property values and some effort must be expended 
in obtaining the best possible thermal properties. The analysis and 
program will continue to be checked as additional flight and ground 
test data becomes available, to both update the thermal property values 
and eliminate the individual approximations and assumptions used in the 
analysis when possible. 
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APPENDIX A 
SAMPLE PROBLEM 


The following sample problem is shown to indicate the form of the 
data input and the program output. A typical charring material sub- 
jected to a constant heating experienced in arc tunnel is presented. 

A sketch of the model is given below: 


v 


q = 95 Btu/sec-ft 


2 



1.5 in. 

I 


0.1 in. 

t 


The various material properties and dimensions are shown in 
table II, program output. The insulation is assumed to be ablation 
material for this problem. The problem coding sheet and subsequent 
data card listing are shown in table 'I. The initial temperature of 
the structure was assumed uniform and equal to 530° R (70° F). Fig- 
ures 8, 9, and 10 are the output data obtained from the plot routines. 
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APPENDIX B 


PROGRAM IN FORTRAN STATEMENTS 




ooooooono 


PRECEDING PAGE BLANK NOT FILMED. 4g 

5IBFTC MAIN A0000 

Aooin 

STRUCTURES and MECHANICS niVlSION A0020 

THERMO-STRUCTURES branch A0030 

THERMAL PROTECTION SYSTEMS SECTION AflOttO 

A0050 

THIS PROGRAM DETERMINES THE PERFORMANCE OF A CHARRING ablator AOOfiO 

A0070 

ANALYSIS AND PROGRAM DEVELOPED Ry DONALD M, CURRY * ES32 AOORO 

A0090 

DIMENSION ESAVEl(S) »ESAVF2<3) »ESAVE3<3) AOlOO 

DIMENSION TITLE (12) rHFADNG( 12 ) »XlDNT < 12» 1 2 ) *TKC<20) »XKC(20) » A0110 

1CPC(20) »TKV(20) ,XXV(20) »TCPV(20),CPV(20) , TIME (30n)#QCON( 300) . A0120 

20RAD(300) »VEL(300) tXNPMd?) .NFPP(12) »NCPR(12) »TXK(20.12),XK(PO*1?) A013O 

3, TCP ( 20 » 12) »CPX(20»12) *RHOBX ( 12) ,XBM( 1 2) ,FMFB(12) *EMBR<12) »HXX< 12) A0140 

4 , GAPX ( 12 ) »FTFST ( 12 ) » BTEST ( 12 ) ,TEmDI (200) ,TX1 (200) *TX2(200) » A0150 

5TX2T ( 10 » 12 ) tTULl ( 200 ) ,TUL2(200> »HX(50) »Tw ( 50 ) » IP ( 50 ) , IR1 ( SO ) , A0160 

6IP2(50) *TUL(50) »IFM(50) * TY (200 ) » A ( 200 ) *B ( 200 ) »C ( 200 ) ,D( 200 ) » A0170 

7R(50) »RHO(50) ,CP(50) ,DXR ( 12) ,XKB( 10*12) »CPB( 10,12) *XMDS(S0) , AOlflO 

BYX(50) *AB(10,l2),RB(10* 12) .CB(10,12) »DB(10»12) *SR(10,12) , A0190 

9RP1 (10*12) ,RB2< 10*12) ,H(12) *S(50) *NPM( 12) A0200 

DIMENSION TTUL(50> *RH0Y1 ( SO ) .RHOy2 ( 50 ) *DrHO<50) ,TCPC(20) A0210 

DIMENSION TlMFP(300) ,PRFS(300) *XC(50) »TX2C(50) #XV(50) ,XDV(50) A0220 

DIMENSION TS(SO) ,SR(S0) A0230 

DIMENSION TTARLE(20).DELTT(20)#IPPC(20) A02U0 

DIMENSION ASAVE1 (3) * ASAVF2(3) *■ ASAVE3 ( 3 ) »BSAVE1 (3) ,BSAVE2<3) , A 0 ? 5 n 

1HSAVE3 ( 3 ) »CSAVE1 ( 3 ) » CRAVE? (3) *CSAVE3(3> *HEAD( 12) * A0260 

1DRAVE1 (3) »DSAVE2(3) »D5AVF3(3) A027O 

DIMENSION XRA ( 30 ) * YA ( 30) A0280 

C A029H 

COMMON TkC * XKC * TCPC * CPC , TK V * XK V * TCPV * CP V , XNPM * RHOBX , XRM , FMBB . A030 0 

1FMFP » NKPB , NCPP , TXK , Xk , TCP ,CPX ,NPM* GAPX tFTFST * RTErT . TEMOI , TX1 , A031 0 

2TX2»TX2T»TUL.TUL1#TUL2»IR»IR1.IR2* A,B*C,D»S*R, AB,RB*CR,DB*SB, A03PO 

3RP1»RR2*TY*RH0Y1»RH0Y2»XMD6*RH0»CP»YK*XKR,CPB#DXB»DT,XL0ST. A0330 

4TABL#TCHAR»TRFCrRHOV,RHOC»FBLDWtFMV#EMC*H300*NKC,NCPC»NKV*NCPV» A0340 

5NP , NMR * NPBS * NPF , TESTg , TFMPI * TXO * TFNV , HENv , FFNV * OLOSS , TLIM » TINT A0350 

COMMON Il»I2,l3*I4»l5*I6.0lN,!NT,DX,XMT*TL»VL,BL,0MP,FRRl,ERR2* A0S60 

1 FPR3 * FRR4 * HV , VPT , CHARK . CHARC » ABLK » ABLC * XmOC » H A0370 

C A03RO 

3000 FORMAT { 12A6 ) A0390 

3001 F0RMAT(1X*12AA> A0400 

3002 FORMAT (6£12»8 ) A0410 

3003 FORMAT (615) A04?0 

3004 FORMAT(115) A0430 

3005 FORMAT(2I5) A04U0 

3007 F0RMAT(2I5*1E14.8) A04R0 

3008 FORMAT (///1X»1 2A6) A0460 

3009 FORMAT (1H1*1X*12A6) A0470 

3010 FORMAT (4E12»8) A0480 

3011 FORMAT(2El2.fl,I6, I5,F13.R*E12.8) A0490 

3012 FORMAT (2£12«8* 16) A0500 

DATA PRV0US/0S4545454B454/ A0510 

RFWIND 11 A0520 

RT0P=9999. A0530 

READ(S*3003)NCASE A0540 

LPL0T=0 A0550 

JCNT=0 A0560^ 
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50 NK = 1 *0570 

I 1 =? *0550 

12=2 *0590 

T?=2 *0600 

14=2 *0610 

T 5=2 *0620 

X6=2 *0630 

T 1 7=2 *0640 

TNT=1 *0650 

XlOST=0.0 *0660 

xMT=0.0 *0670 

XMDT=0.0 *0650 

FRR 1=0 , 0 *0690 

F'PR2=0 • 0 *0700 

FPR3=0 , 0 *0710 

FRR4=0,0 *0720 

ICT=0 *0730 

TCONTrO *0740 

XMDC=0.0 *0750 

NXP=1 *0760 

XlSTV=0,t) *0770 

MPS=2 *0750 

FPR5=0 » 0 *0790 

TPCT=0 A0500 

TTTP=0 A051O 

TPL0T=1 A0520 

NXA=1 *0530 

NXB=1 *0540 

NXC=1 *0850 

NXD=1 *0560 

S*VY3=-ll>0, *0870 

5*VY4X=-100. A0880 

5X0=0, 0 *0890 

500T=0,0 *0900 

*0910 

GFNERAL llTl.F OF PROflLEM *0920 

100 rFAD(5*3000) (HEAD(K) ,K=1 *12) A0930 

WRITE (6 ► 3009 ) ( HEAD <K ) , K=1 . 12 > *0940 

LPLOT=LPLOT+l A09RO 

WRITE (11) (HEAD(I) * 1=1* 12) *0960 

wP ITE (6*110) *0970 

110 F0RMAT(//1X,11HINPUT DATA,//) A0950 

RF AD ( 5 , 3000 ) (TITLE(l) *L=1 *12) *0990 

IF (TITLE ( 1 ) »F0 , PRVOUp) GO TO 150 *1000 

HFAD ( 5* 30.1 1 ) TL, IM, TINT* NPTT , NPLOT , OMP»TDMp *1010 

RFAD(5*3ul2) (TTARLF(T),0FLTT(I),IPRC(I),I=1.NPTT) *1020 

T=TINT *1030 

nTS=DFLTT(l) *1040 

DT=DELTT ( 1 ) /3600 . 0 *1050 

WRITE (6* 120 ) TLIM,TINT*NPTT *1060 

120 FORMAT(1HO*11HTIMF LIMITS, 1PF10. 4, 4X,I3HTNITI*L TIMF=,1PF10.4,4X.5 *1070 

1HNPTT=,I4) *1080 

WRITE (6* 122) *1090 

122 FORMAT (//5X,4MTIMF*10X,9MTIMF STfP*6X* 13HPRINT CONTROL ) *1100 

WRITE (6* 124 ) ( TTAPLE ( I ) , DFLTT ( I ) , I PRC (I), 1 = 1* NPTT ) *1110 

124 FORMAT (5X * 1PE1 0,4* 6X , 1 PF 10.4 »OX , 14 ) *1120 

C *1130 
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c location factors for convfctive and RAOIATIVE HEATING A1140 

ISO RFAD(B,3000) (TITLE(L) »L=1 .12) *1150 

IF (TITLE ( 1 ) #EO»PRVOUSl GO TO 200 A1160 

RFAD(5,3002) FCONV.FRAD A1170 

WRITE ( 6 * 1S5) FCONV.FRAD A1180 

1S5 FORMAT ( 1H0 .6HFC0NVZ. 1PE12 ,5»4X5HFRADs» 1PF12.5/ ) A1190 

A1200 

PROPERTIES OF ABLATION MATERIAL A1210 

200 RFAD(5,3000) (TITLE(L) ,L=1»12) *1220 

TF(TITLE(1) .EO.PRVOUS) GO TO 300 A1230 

RFAD(5»3000) tHFADNG(K) ,K=1,12) A12aO 

RFADt 5,3002) TABL,TCHAR»TREC,PHOV»RHOC»FRLOW,FMV,FMC,H30n»VL,HV, A1250 

1 VPT » FV , T V * CHARK » CHARC * ABLK * ABLC A 1260 

RF A ( 5, 3003) NP»NKC »NCPC , NKV »NCPV» NREC A1270 

READ ( S , 3o02 ) (TKC(K) ,XKC(K) .Ksl.NKC) A1280 

RFAD<5,3002) (TCPC(M) ,CPC(M) ,M=1,NCPC) *1290 

RFAD(5,3002) (TKV(L) ,XKV<I ) »L=1»NKV) A1300 

RFAD(5,3002> (TCPV(N) ,CPV(N) ,Nsl,NCPV) A 1310 

RFAD(5»3U02) asm ,RR(I) »I=1,NRFC) A1320 

WRITE ( 6 , 3008 ) < HFADNfi ( K ) , KZl , 12 ) A1330 

WRITE (6,210) TABL»TChAR»TREC»RHOV*RHOC*FRLOW,FMV,EMC,H300»VL.HV. A13U0 

1vPT,FV,TV,CHAPK,CHARc»ARLK»ARLC A13S0 

210 FORMAT (1H0»5HTABL=» 1PF12.5»3X,6HTCHAR=, 1PF12.5»3X»5HTREC=, 1PF12.5, A1360 

13X,5HRH0v=,1PF12.5,3x,5HRH0C=,1PF12.S,21X/1X,6HFbL0W=,1PF12.5,4X,4 A1370 

2hFMV=,1PE12.R.4X,4HEmC=.1PE12.5*SX*5HH300S*1PF12.5.5X.3HVL=»1PE1R. A1380 

35/4X » SHHV=, 1PF12.5.4X, 4HVPTS, 1 PFl 2 ,5, 5X , 3HFVs* 1 Pf 12»5,5X» SHTVr * 1PE *1390 

llP.5,2X,6HCHARK=,lPFl2.R/tX»6HCHARC=,lPEl2.5,3X»RHARLK=,lPE12.5,3X A1400 

2,5HABLC=,lPFl2.5/) A1410 

VLISVL A1420 

VL=VL/12.0 A1430 

VLV=VL A1440 

WPITE(6,2P0) NP » NK C » NCPC • NKV • NCPV , NREC A14R0 

220 FORMAT ( 2X , 3HNP= , 1 14 » aX , 4HMKC= , 1 la , 4X . 5HNCPC= , 1 ia , 4X , aHNK V= » 1 la » 4X , A1460 

1 RHNCPVS » 1 14 , ax » 5HNRECS » 1 14 ) A 1470 

WRITE (6 * 221 ) A1480 

221 F0RMAT(/32X»1SHVIRGIN MATFRI AL/20X . 7MTHERMAL. 38X , 8HSPFCIFIC/3X » 1 1 H A1490 

1TFMPERATURE,4X» 12HCOnDUCTTVITY» 19X* 11HTEmPERATURF»7X,4HHfAT) A1500 

KLLL=MIN0 (NKV. NCPV) *1510 

WPITE<6»222) (TKV(L) ,XKV(L) »TCPV(L) »CPV(L) »L=1 .KLLL) A1520 

222 f 0RMAT(2X»1PF12.5»4X, 1 PE12.5 . 1 8X , 1PF12 .5 .3X,lPFl2,5) A1R30 

IF(NKV-NCPV) 223.227,225 A1540 

223 KLLLL=KLLL+1 A1550 

WRITE ( 6 , 224 ) ( TCPV(L) ,CPV (L) .LsKlLLL.MCPv) *1560 

224 FORMAT ( 48X , 1PF12 . 5 » 3x » 1PF1 2.5) A1570 

GO TO 227 A15B0 

225 kLLLL=KLLL+1 A1590 

WRITE (6 .226 ) ( TK V (L ) , XKV (L ) »L=KLl.UL»NKV ) A1600 

226 FORMAT (2X»1PF12#5,4X#1PE12,5) A1610 

227 WRITE(6.2P8) A1620 

228 FORMAT (//33X.14HCHAR MATFRI AL/20X.7HTHERMAL. 38X, 8HSPFC IF IC/3X. 11 H A1630 

lTFMPERATURE»aX.12HC0NDUCTIVITY,l9X.llHTFMPERATURF,7X,aHHEAT) A1640 

KLLL=MINO(NKC,NCPC) A1650 

WRITE(6,222) ( TKC <L > , XKC <L) .TCPC(L) #CPC <L) ,L=1 ,KlLL) A1660 

IF(NKC-NCPC) 230,235,232 A1670 

230 KLLLL=KLLL+1 A1680 

WRITE(6»2P4) (TCPC(L) ,CPC(L).L=KLLLL,NCPC) A1690 

GO TO 235 A1700 



on no 


52 


232 KLLLLsKLLL+1 *1710 

wRITE(6.226) <TKCtL) , XKC (L> .LsKLlLL.NKC) *1720 

235 WRITE(6»240) *1730 

240 FORMAT <//28X,?3HsuRFaCE RFCESDION TABLE//25X, 1 IHTFMPERATijRE. 8X» 11 H *1740 

1DR - TN/SFC) A1750 

WRITE (6*245) (TS<I),SR(I>,l31,NRFC> A1760 

245 FORMAT (24X. 1PF12.5»7X»1PF12.5) A1770 

*1780 

PROPERTIED OF TRAJECTORY A1790 

300 RFAD(5,3000) (TITLE<L) »Ls1,12> *1200 

IF (TITLE (D.FO.PRVOUO 60 TO 400 *1210 

RF AD ( B» 3000 ) ( HEADN6 ( L ) » Lsl » 12 ) *1220 

RFAOtD.3004) NTRAPT *1830 

READIES!)! 0) (TlMF(K) ,QCON(K) ,QR*D(K) ,VEL(K) ,KSl, NTRAPT) A1840 

WRITE (6. 3008) <HE*DNfi< L) ,L=1 . 12 ) A1850 

wRITE(6»310> NTRAPT *1260 

310 FORMAT ( 1H0 » 27H NO. OF TRAJECTORY POINTS 5,114) *1870 

wPITE(6»3?0) A1880 

320 F0RMAT(//PX»4MTIMF»AX»12H0 C0NVEcTIVE.4X.I1H0 RADIATIVE. 7X.RHVELOC *1290 

1 ITY ) *1900 

WRITE(6»330> ( TIME (K ) . QCON(K ) .OR*D(K) . VEL(K) .KSl, NTRAPT) *1910 

330 FORMAT ( 1R4E16.5) *1920 

*1930 

PROPERTIED OF BACK-UP STRUCTURE *1940 

400 RFAD(D.3000) (TITLE(L) .Lsl ,12) *1950 

IF (TITLE ( J ) .EO.PRVOUO 60 TO BOO *1960 

WRITE(6»4I0) *1970 

410 FORMAT (//10X.31H PROPFRTIFS OF BACKUP STRUCTURE/) *1980 

RFAD(B»3007) NMP.NPBs.BL *1990 

RFAD ( R» 3002 ) ( XNPM ( K ) . Ks 1 . NMB ) *2000 

RFAD(5,41D> ( NKPB ( I ) , NCPB < I ) . Isl ,NMB) *2010 

415 FORMAT(10I5> *2020 

no 420 Ksl.NMR *2030 

NPM(K)=XNPM(K)+0. 00000002 *2040 

420 CONTINUE A2050 

WRITE (6. 425) NMP.NPBS.BL *2060 

425 FORMAT (/4X.3SHN0. OF MATFRIALD IN BACK-UP SHIFLOs, 1I4/4X.40HT0TAI. *2070 

INUMBER OF NODFS IN RaCK-IIP SHIELDS, 1I4/4X.28HTHICKNESS OF BACK-UP *2020 

29HIFLD=,1PE1?.5//) *2090 

PLsBL/12 • 0 *2100 

no 440 Isl ,NMP *2110 

LKsNKPB(I) *2120 

LCPsNCPB(I) *2130 

RFAD ( S, 3000 ) ( ( XIONT ( K , I ) ) »Ksl . IP ) *2140 

rFAD(B,3o02) ( <TXK(J,I) ,XK(J,T) ) . J=1,LK) A2150 

RFAD(D,3002) ((TCP(J,I),CPX(J,I)),Js1 # LCP) *2160 

wPITE(6»432) (XIDNT(K.I) .KS1.12) *2170 

432 FORMAT ( //I 2*6 ) *2180 

WPITE(6»433) *2190 

433 FORMAT (//20X,7HTHFRMaL»38X.8HDPEcIFIC/3X,11HTFMPFRATURE.4X»1PHCOND *2200 

lltCTlVlTY » 19X, 1 1 HTFMPfPATuRE,7X»4HHE*T) * 2210 

KLLLsMIN0(LK,1 CP) *2220 

no 434 Nsl.KU L A2230 

WRITE(6,222) ( TXK ( N, T ) , XK <N, I ) .TcP(N, I ) ,cPX (N, I) > *2240 

434 CONTINUE *2250 

IF(LK-LCP) 432,440.437 *2260 

435 kLLLL=KLLL+1 *2270 
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no 43ft N=KLLLL»LCP 

WRITE (6# 224) <TCP<N»T> ,CPX(N,T) ) 

436 CONTINUE 
GO TO 440 

437 kLLLLsKLLL+1 

no 43ft NsKLLLL »lk 
WRITE<6»226) (TXK(N,T) »XK(N.I) ) 

438 CONTINUE 
440 CONTINUE 

RFAD(R»3il02) (RHOPX(L) »XBM(L) #EMFB(L> ,EMrR(L) »L=1 »NMR) 

RFAD(5#3o02) <H(J) ,RAPX(J) .FTERT(J) #BTFST< J) . Js1#NMR) 

wPITE<6#4«50) 

4*50 FORMAT <///55x.10HEMI«;*;IVTTY/8X.RHMATEPIAL»5X.7HDeNSITY# 7X#9HTHICKN 
lFSS » 7X » 5HFR0NT » 9x ► 4HB ACK , 7X » 14HNODES/MATFR I AL/ ) 
no 460 LLJSl.NMP 

WRITE ( 6 ► 4R5 ) LLJ.RHOrx(LLJ) #XPM<LLJ> »FMFR(LLJ) »EmRB(LLJ) .XNPM(LL j) 
4*55 FORMAT(UX#lI1 # 8X # 1PF1 0 . 4 #4X , IPEl 0 .4 »4X# 1 PE10 .4 »4X » 1PF10.4 #6X . 1PF1 
10. 4/) 

4ft0 CONTINUF 

WRITE(6#465) 

4ft5 FORMAT ( //4X.ftnHADOITtONA| DATA FOR INDIVIDUAL MATERIALS IN BACKUP 
1FTRUCTURE//11X»RHMATFRIA| ,5X.!6HFILM COEFFICIENT » 5X» 13H6AP THTCKNF 
2sS , 8X , 5HFTEST . 13X » 5HRTE5T ) 
no 48n j=i#nmr 

wRITE<6»470) J# H(J) .GAPXCJ) #FTFST< J) #BTEST(J) 

470 FORMAT <13X#1 13# 1 2X # lPElfl . 4»9X . 1PF10 .4 »7X, 1PE11 .4,7X»1PE1 1.4/) 

4R0 CONTINUE 

PROPERTIES OF ENVIRONMENT 
500 RFAD(*5#3000> ( TITLE (L) * Lsl . 12 ) 

TF (TITLE ( 1 ) .FO.PRVOIK) 60 TO 600 
RFAD(R#3o00) (HEADNG(L) »L=1»1?) 

RFAD(«5»3002> TENV»HEnV#FFNV»QL 059 
wRITE<6#3008> (HFADNs(L) .L=l »12) 
wRITE(6»b?0) TENV»HEnV*FFNV»QLOS«; 

520 FORMAT </4X>l2HTFMPERATUPF=»lPF12.5»4X,17HFILM COEFFICIENTS# 1PF12.S 
1 .4X» 12HV1FW FACTORS# 1PE12 ,5»4X»7H0 LOSTs, 1PE12.5) 

INITIAL TEMPERATURE DISTRIBUTION 
600 PFAD(S#3000) ( TITLE (L )# Lsl . 12 ) 

TF(TITLE(1) .FO.PRVOUF) 60 TO 700 
RFAD05.3000) (HEADNG(L) #Lsl#lp) 

npfsnp+npps 

TLsVL+BL 

XNPSNP 

nxsVL/(XNP-i.n) 

nXXsDX 

RFAD(S#3002) TEST2»TFMPI,TX0 
IF ( TE5T2 ) 6l0»620»62n 
610 PFAD ( *5# 3002 ) (TEMDI <<) #K = 1 »NPF) 
no 61*5 Ksl »NPF 
TX1 (K ) STEMDI (K ) 

TX2(K)STX1 <K> 

TULKK)STXl(K) 

TUL2(K)sTX1(K) 

615 CONTINUE 
I.sNP+1 


A2280 
A229 ft 
A2300 
A2S10 
A2320 
A2330 
A2340 


A2350 

A2360 

A2370 

A2380 

A2390 

A2400 


A24t0 


A2420 
A2430 
A2440 
A2450 
A2460 
A2470 
A2480 
A24Q0 
A2500 
A2510 
A2520 
A2530 
A2R40 
42550 
A2560 
A2570 
A2580 
A2B90 
A2600 
A2610 
A2620 
A2630 
A2640 
A2650 
A2660 
A2670 
A2680 
A2690 
A2700 
A2710 
A 2720 
A2730 
A2740 
A2750 
A2760 
A2770 
A2780 
A2790 
A2800 
A2810 
A2820 
A2830 
A2R40 
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no 619 1=1 » NmP 




A2850 


) N=NPM(I ) 




A2860 


no 617 J=1 »LN 




#2870 


T*2T< J,I )=TFMni(L> 




A2880 


1 =L+1 




A2890 

617 

CONTINUE 




A2900 

bl 9 

CONTINUE 




A2910 


00 TO 625 




A2920 

620 

CALL TEMPO 




A2930 

625 

wPITE(6»3008) (HEADNg(L) ,1=1.12) 




A2940 


TF(TEST2) 630, 635,640 




A2950 

630 

WRITE<6»632) 




A2960 

632 

FORMAT (4X *52HTEMPERATllRF DISTRIBUTION IN HEAT 

SHIFLD 

TS 

apbritary/ 

A2970 

1) 




A2980 


wPITE (6,633) ( TFMOI (K) » K=1 ,NPE ) 




A2990 

b33 

FORMAT(lP#E12,5) 




A3000 


C,0 TO 64 b 




A3010 

635 

wPITE (6.637 ) TEMPI 




A3020 

637 

FORMAT C//4X.64HTEMPFRATUPF DISTRIBUTION TN HEAT SHIFLD 

IS UNIFORM 

A3030 

IaND EOUAL TO ,1PE10.u/) 




A3040 


00 TO 64b 




A3050 

6U0 

wPITE (6.641) 




#3060 

toUi 

FORMAT (4X.54HI INEAR TEMPERATURE DISTRIBUTION 

ASSUMED 

IN 

hfat chiel 

A3070 

m/> 




A3080 


wRlTE(6»h33) (TEMOI(L) .Lsl.NPE) 




A3090 

6b5 

TF(DMP) 700,700,646 




A3100 

bub 

WRITE (6 , 647 ) 




A31 10 

647 

FORMAT (//> 




A31?0 

648 

wPITE (6, 649) (TXKL) ,TX2(L > »L=1»NPF) 




#3130 

6U9 

FORMAT ( 2X , 1PE1 2*5. 4X ,1 PE 12,5) 




A3140 


WRITE(6»650) 




A3150 

650 

FORMAT (//) 




A3160 






#3170 


enthalpy AS A FUNCTION OF TFMPERATURE 




A3180 

700 

READ (5,3000) (TITLE(L) »L=1 ,12) 




A3190 


TE(TITLE(1),E0.PRV0U5) 00 TO 725 




A3200 


«EAD(5,3004> NHP 




A3210 


RFAD(5,3002) (HX(K) ,TW(K) ,K=1,NHP) 




#3220 

7?5 

00 720 1=1 , Np 




#3230 


TP ( I > =0 




#3240 


T P 1 ( 1 ) =0 




A3250 


TP?<I)=0 




A3260 


TEM( I )=0 




A3270 


XMDO ( T ) =0 , 0 




A32B0 

728 

CONTINUE 




• 3290 


wPITE (6,730 ) 




#3300 

730 

FORMAT ( 1H1 »12H0UTPUT DATA.//) 




#3310 


XC(1)=0.0 




A3320 


no 740 I=2»NP 




#3330 


XC(I)=XC(I-1)+DX 




#3340 

7U0 

CONTINUE 




#3350 

730 

IF(T-TIME(NK) ) 765,770,760 




A3360 

760 

nK=NK+1 




#3370 


TF (NK-NTKAPT) 750,750,762 




#3380 

762 

WPITE(6,763) NK 




#3390 

763 

format (1H0,33H THE VALUF OF NK IS IN ERROR, NK=»1I4) 



#3400 


00 TO 905 




#3410 



on o o no no 
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765 TF(NK-2) 762,766,766 A3420 

766 OC0NX=QC0N(N«-l)+( ( OcON< NK >-QCON < NK-1 ) > / ( TIME ( NX ) -TIMF < Nk-1 ) ) ) A3430 

1*(T-TTMEINK-1) ) A3440 

OC0NX=FC0NV*QC0Nx A3450 

0RADX=QRAD(Nk- 1)+( (ORAD(NX)-ORAD(NK-l) )/(TIME(NK)-TlMF(NK-l) ) ) A 3460 

1 * ( T— T T MF ( NK- 1 ) ) A 3470 

ORADX=FRAD*GRADX A3480 

VFLX=VEL(NK-1 )+< (VEL(NK)-VEL(NK-I ) ) / ( TIMf ( NK )-TImF ( NK-1 ) ) ) AS4Q0 

i*(t-ttmf(nk-h ) assoo 

60 TO 775 A3510 

770 0C0NX=FC0NV*0C0N(NK) A3520 

ORADX=FRAD*GRAD(NK ) A3530 

VFLX=VEL(NK) A3540 

A3550 

COMPUTE HFAT RLOCKAGf AT FRONT SURFACE A3560 

775 TF( 117-1) 778,776,776 A3570 

776 TF ( I 17— NHP ) 777,777,778 A3580 

777 TF(TX?(lNT)-TWUl7) ) 782,788,780 A3590 

778 WRITE ( 6 » 779 ) TXP(INT) A3600 

779 FORMAT <1HO»80H THF RANGF OF THE FNTHALPY-TEMPFRATURE CURVF FIT WAS A3610 

1FXCFEDED AT A TEMPERATURF 0F,1E10,4) A3620 

60 TO 905 A 3630 

780 T 1 7=1 17+1 A 3640 

60 TO 776 A 3650 

782 TF(TX2(INT>-TW(I17-1) ) 784,788,786 A3660 

784 T 17=1 17-1 A 3670 

60 TO 775 A3680 

786 HWsHX ( 117-1 ) + ( (HX(I17)-HX (117-1) )/(TW(Il7)-TW( 117-1) ) ) A3690 

1*(TX2( INT )-Tw( 117-1) ) A3700 

60 TO 789 A3710 

788 HW=HX(I17) A3720 

789 mTX=H300+( (VFI X**2 > /60056 .5 ) A3730 

ORL0CX=(FRL0w*XMDGClNT)*(HTX-HW) )/3600,0 A3740 

A3750 

COMPUTE HFAT TN DUE TO SURFACE COMBUSTION A3760 

XMD0=XMDC A3770 

CALL OXIDAT(XmDO,OOXtO) A3780 

A3790 

COMPUTE Q-HOT WALL A3800 

TF(TDMP,EO.O. ) GO TO 4001 A3810 

TF(T.GE.TDMP) DMprl.n A3820 

4001 7=(HTV— Hw)/ (HTX— H300) A3830 

IF(Z-I.O) 790,792,793 A3840 

790 TF(Z) 791,791,793 A3850 

791 oHw=0 . 0 A3860 

60 TO 1790 A3870 

792 OHW=QC0NX A3880 

60 TO 1790 A3890 

793 OHw=Z*QCONX A3900 

1790 7ZZ=(0HW-0BL0CK)/0HW A3910 

TF(ZZZ-0.2) 1798,1798,1794 A3920 

1708 OPLOCK=0,8*QHW A3930 

A39U0 

NFT HFAT INTO FRONT SURFACE A3950 

1794 IF(IEMCINT)) 795,795,797 A3960 

795 TF(TX2(INT)-TCHAR) 746,796,797 A3970 

796 FMX=EMV A3980 
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GO TO 79« A3990 

707 1FM(INT)=1 A 40 oo 

FMX=EMC A4010 

708 OlN=0RADX+QHw+Q0X ID-fjRLOCK- ( 4 .83.33E-13 ) *FMX*FV* ( { TX2 ( TNT ) **4 ) - A4020 

1 ( TV**4 ) ) A4030 

if(Dmp) ao4»ao4»aoo A 4040 

800 wPITE(6i80l> A4050 

801 FORMAT! ///) A4080 

WRITE (6»an2) OC0NX,ORADX,VELX,HTx,HW.7»0RL0CK,QHw#Q0XTD,GTN A4070 

802 FORMAT!) X»6Hqc 0NX=, 1PFl2.fi. 2X,6HQPADX=»1PF12.S»2X.5HVFLX=»1PF12. *5# A4080 

l7X,4HHTX=,lPei2.5,2X,3HHw=»lPF12.5/lX,2H7=»lPF12.5»2X,7H0RL0CK=,lP A4000 

2F12.b,2X»AH©HW=.lPEl2.5»7X,6HOOXlD=»lPEl 2. 5. 2X.4H0INs.lPF12.fi/) A4100 

804 OTN=QTN*3800 . 0 A4110 

A4120 

CHECK FOK FROMT SURFACE RFCESSION (CHAR I.AYFR REMOVAL) A4130 

CALL RECESSCXmDC.XLOsT.TRFC.DT.RHOC.TS.SR. TX 2(1) .NREC.NRS.ERRS.SXO A4140 
l.SDOT.DMP) A4150 

tf(Err 5) 80fin.80so.9ns A4ian 

BOfiO VLV=VLV-XLOST A4170 

XI STV=XLSTV+Xl OST A4180 

XI STI=XLSTV*12.0 A4190 

nXV=VLV/ (XNP-1 .0) A4200 

XV(1)=0.0 A4210 

no 17A0 I=2»NP A4220 

XV (I )=XV(T-1)+0XV A4230 

17A0 CONTINUE A4240 

nXsOXV A4250 

IF (ERR4 ) 806,808, 80S A4260 

80b GO TO 90b A4270 

806 CALL COFFFINPFT.snOT) A4280 

IFIDMP) 8069, 8069, 80ftl A4290 

8081 WRITE 1 6» 8062 ) A4300 

8082 FORMAT ( /IX » 23H COEFFICIENTS FOR SWUFT/) A4310 

00 8086 I=1»NPFT A4320 

WPITE (6.8064 ) A(I), RII) ,C(I),D(I),I A4330 

8084 FORMAT !1H0.5ha(I)=» 1PF12.S.2X ,5HR( I ) = * 1PF12 .5 » 2X, 5HC (T)=»1PE12.S.2 A4340 

lX .5HD (I)=»1 Pf12.5»2X.7HI=*I3) A43S0 

8086 CONTINUE A4360 

8089 TFIFRR2) 807,807,80*5 A4370 

807 TFIERR3) 810,810,808 A4380 

608 wPITE(6*809) TKK A4390 

809 FORMAT ( 1H0 » 18M ThF VaLUF OF IKKS.1I4) A4400 

GO TO 90b A4410 

810 CALL SWUFT(A,R.C,D»Ty»NPFT.DMP) A4420 

827 no 828 1=1 .Np A4430 

TX1 ( I ) =TX2 ( I ) A44UO 

TX2 ( 1 ) =TY ( I ) A4450 

828 CONTINUE A4480 

CALL OON2 ( XLOfiT » XV »Tx2» NP .XC *TX2C* XDV.KKv. XLSTV * OXX) A4470 

830 CALL ABLATE A4480 

XMDT=XMOG(INT)+XMnC A4490 

lT=NP+l A4fi00 

00 1815 1=1 ,NMB A4510 

1 LT=NPM ( 1 ) A4520 

IF(I.FQ.I) 60 TO 1812 A4530 

TFIGAPXI 1-1) .FO.O, ) GO TO 1812 A4540 

KKT=1 A4B50 
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GO TO IBIS 

1812 KKT=2 

1813 00 1815 JrKKT.LLT 
TX2T(J»I)=TY<l.T) 

I T=LT+1 

1815 CONTINUE 

no 1819 I=1*NMB 
tF(I.FQ.I) GO TO 181ft 
TF(GAPX(I-1).FQ.0.) GO TO 1817 
GO TO 181«> 

1816 TX2T ( 1 ,1 )=TY(NP) 

GO TO 181Q 

1817 l. XsNPM (1 — 1 ) 

TX2T(1,I)=TX2T(LX,I-1) 

1819 rONTINUF 
I PeNP+1 

no 833 1=1 * NMP 

t. 7=NPM(I ) 
nO 833 J=1 *L7 
TX2(LM>=TX2T(J»I) 
i M=LM+l 
833 CONTINUE 

nO 5834 I=2*NPTT 

TF (T-TTABLE( I ) ) 5835,5835*5834 

5835 nTS=DFLTT(I-n 
TPRCT=IPKC(I-1 > 

OTsOELTT (t“l)/3600»0 
GO TO 5838 

5834 CONTINUE 

nTS=DFLTl (NPTT) 

TPRCTsIPKC(NPTT) 

OTrDELTT (NPTT 1 /3ft00 . n 

5836 TCT=ICT+1 
5838 VLTEM=SAVY3 

CALL TS0THM(XV*TX2*l060. *NP»SAVEIT> 

GAVEIT=SAVEIT+XLSTV 

TF(SAVY3.LT,SAVEIT)SAVY3=5AVEIT 

IF(VLTEM.FQ,SAVY3)G0 TO 839 

SAVX=T 

GAVYIsXLSTI 

SAVY2=TX2(NP> 

CALL ISOTHM ( XV * TX2 * 1460, ,MP*SAVY4> 

839 Bl. TFM=SAVY4X 

CALL TSOTHM(XV»TX2>l460 . *NP,WFKFfP) 

wFkEEP=WEKEEP+XLSTV 

T F ( SAVY4X.LT . WEKEFP ) SA VYuXsWEKEFP 

TF(BLTEM.FQ.SAVY4X)G0 TO 838 

SAVEXX=T 

SAVY1X=XLSTI 

SAVY2X=TX?(NP1 

CALL TSOTHM( XV* TX2 * 1 OftO. *NP»SAVY3X) 

838 CONTINUE 

IF(IPPCT-ICT) 835*835*840 

835 WPITE(6*837> T * QCONX , GRAnx * VFL X * XMDG (INTI* XPOC . X mOT * XLST I » QHW 
837 FORMAT(1HO*5HTIME=* 

1 1PE12.5.2X,12H0C0NVFCTTVE=.1PF12.5,?X»11H0RA0IAT 

1TVE=*1PF12.5*5>X*9HVFlOCITY=*1PE12.5/1X*1AHGAS APlATION RATE=,lPFl? 


A4560 

A4570 

A4580 

A4590 

A 4600 

A4610 

A4620 

A 4630 

A4640 

A4650 

A4660 

A4670 

A4680 

A4690 

A4700 

A4710 

A4720 

A4730 

A4740 

A4750 

A4760 

A4770 

A4780 

A4790 

A4800 

A4810 

A4820 

A4830 

64840 

A4850 

A4860 ! 

64870 ! 

A4880 

A4890 ; 

64900 i 

A4910 | 

A4920 

A4930 ; 

64940 I 

A 4950 I 

6496 0 ; 

A4970 I 

A4980 

A4990 f 

65000 t 

65010 i' 

A5020 

A5030 ; 

A5040 : 

A5050 '{ 

A5060 V- 

A5070 ; 

A5080 ?: 

A5090 fe 

A5100 U 

A5110 ! 

65120 •; 
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2.S.2X.1QHCHAR ablation PATEs,1PFi2.5»?X.?OHTOTAL ablation RATES . 1 P A51S0 
3F)2,fi/lX.16HRFCFSSI0N DFPTHs, 1PF)?.5.?X» 1 OHGHOT wALLr, lPF12.fi) A5140 

840 TST+DTS A5150 

841 tF (NPI..OT .ME, t ) 00 TO A4? A5160 

CALL SAVE ( ASAVE1 . ASAv/F2» ASAVF3»UfiFA* NXA » XL STI * DTfi » TLIW.T. VALllFA ) A5170 

CALL fiAVE<BSAVEl,RSAvF2,RSAVF3.UfiFB.NXB»TX2<NP> .DTS.Tl IM.T.VAI UFR ) A5180 

CALL TSOTHM<XV»TX?.l060.,MP,Y3> A5190 

CALL PAVE (CSAVEl ,CSAvF2,CSAVF3.UfiFC.NXC.Y3»DTfi»T| IM.T.VALUEC) A5200 

CALL TSOTHM ( XV. TX? » 1460. *NP » YU ) A5210 

CALL SAVE(DSaVE1 ,DSAvF2.PSAVF3.USFD.NXD»YU,DTS.TlTM»T, VALUED) A52?0 

TFfUSFA.NF. 0.0)00 TO 9842 A5230 

TFtUSFB.NF. 0.0)00 TO 984? A52U0 

TF(USFC.NF.O.O)GO TO 984? A52S0 

TF(USFD.NF. 0.0)00 TO 984? A5260 

O 0 TO 9843 AS?70 

9842 XPL0T=T-0TS A5?80 

YPL01 isVAl UEA A5290 

TF<USFA.NF.0.O)YPLOTl=UfiFA A5300 

YPLOT?sVALUER A5310 

I F ( USFB , NF . 0 . 0 ) YPLOT?=UfiFF A53?0 

YPL0T3SVALUEC A5330 

T F ( USFC . NF . 0 . 0 ) YPL0T3SU5EC A53U0 

YPLOT UsVALUFD A53B0 

T F ( USFD . NF . 0 . 0 ) YPLOT usUSFO A5360 

wPITfc ( 1 1 )XPlOT»YPL0T) » YPI 0T2.YPI OT3.YPLOT4 A5370 

9843 TF(ICTP.NF.O) GO TO 842 A5380 

TCTPs) A5390 

XPLOTsT A5400 

YPLOTIsXLfiTI A5410 

YPLOT?sTX?(NP) A5420 

CALL TSOTHM(XV»TX?»l060. »MP.YPL0T3) A5430 

CALL TSOTHMCXV.TXP. 1460. .MP.YPL0T4) A5440 

wPlTE (11 ) XPloT » YPLOT 1 » YPL0T2 » YPL0T3. YPL0T4 AfiUfiO 

842 TF(IP»CT-TCT) 84&.845.900 A5460 

845 WPITE (6. 8fiO ) T A 5470 

TPCTsTPCT+1 A5480 

TF(IPCT.E0.2)TPCTSO A 5490 

TF ( IPCT • EO • 0 ) TCTPsO A5500 

850 FORMAT (1H0.7PMTFMPERATUPF DISTRIBUTION In HFAT SHIELD AT THF FNP 0 A5510 

IF THE TIwF STFP, Ts , 1 PF) ? .5 » 1 X .7HSEC0NDS// ) A5520 

wP ITE ( 6 » 880 ) A5530 

880 FORMAT (4X.49HTEMPFR A TURF DISTRIBUTION IN THF ABLaTINQ MATERIAL//) A5540 
kK V=KKV+ 1 A55fi0 

wP ITE ( 6 » 882 ) (TX?C(I) ,ls) ,KKV) AfififiO 

882 F0RMAT(6X. 1PF12.5. 1P5F16.5) A5570 

TJsNP+1 A5580 

wP ITE (6. 864 ) A5590 

864 FORMAT (//4X.40HTRMPFPATURF DISTRIBUTION TN THF PACK-IJP STRUCTURF// A5600 

1) A5610 

WPITE(6»862) ( TX2 ( I ) » ISI J.NPF ) A5620 

wPlTE ( 6» 865 ) A5630 

865 FORMAT (//) A5640 

TCT=0 A5650 

900 CONTINUE A 5660 

TF(T-TLIM) 7fiO.750.9nfi A5670 

905 tF(NPLOT.NE.I) GO TO 909 A5680 

XAVY3=SAVY3-SAVY1/12, A5690 
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XAVY4X=SAVY4X-SAVY1X/12. A5700 

T F ( S AVX • EG » XPt OT ) GO TO 9n05 A5710 

wPITE(ll)SAVX,SAVYl»SAVY2.XAVY3»SAVY4 A5720 

9005 TF<SAVEXX.EG.XPLOT)GO TO 9006 A5730 

«;AV41=SAVY4X*12. A5760 

9006 FAV3I=SAVY3*12. A5750 

WRITE ( 1 1 ) SAVEVX*SAVY 1 X» SAVY2X #SAVY3X > XAVY4X A57U0 

WRITE (6»929)SAV3I rSAv4l A5770 

929 FORMAT ( 1H0 » 23HMAXTMUM 1060 ISOTHERM =E16.8,2X23HmAXIMUM 1460 ISOTH A5780 

lFRM =F16.A> A5790 

WRITE 1 11 ) STOP » ST OP » STOP » STOP* STOP A5800 

909 TR (LPLOT .NE.NCASEIGO TO 91 1 A5810 

DATA FND/6H FND / A5820 

wP I TE ( 1 1 ) END , FND . END , FNn , END » FND , END » END . FND . FND . FND , FND A5830 

OUIT=A8A». A5840 

WRITE ( 11 )OUIT,OUIT#Ol)IT»OUIT,OUlT A5B50 

FMD FILE 11 A 5860 

RFwIND H A5870 

911 TF(TESTP) 9l0*930»93n A5880 

910 no 920 JJK=1,MPF A5890' 

TXKJJK>=TEMDT(JJK) A5900 1 

TV2(JJK)=TX1(JJK) A5910 

TUL1(K)=TX1<K) A5920 

TUL2(K)=TX1 <K > A5930 

920 CONTINUE A5940 

TL=NP+1 A5950 

no 926 I=1#NMP A5960 

Tl NrNPM(l) A5970 

DO 924 JrltlLN A59S0 

TX?T(J»I)=TEMni(IL> A5990 

TL=IL+1 A6000 

924 CONTINUE A6010 

926 CONTINUE A6020 

fiO TO 940 A6030 

930 CALL TEMPO A6040 

940 T=TINT A6050 

nX=DXX A6060 

nTS=DELTT(l) A6070 

DT=DELTT(1 )/.3600.0 A6080 

VI V=VL A6090 

fiO TO 50 A6100 

END A6110 
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SIBFTC (“OFF 

C THIS SUBROUTINE DFTERMINFS THF COEFFICIENTS OF THE MATRIX 

SUBROUT INF COFFF ( NPFt , SHOT ) 

C 

DIMENSION TITLE ( 1? ) »HEADNG( 12) t X IDNT (12,12), TKC(20)»XKC(?0> , 
lCPC<20) ,T*V(?n) ,XKV(? 0 ) ,TOPV(?0) ,CPV(20) ,TIME(30n) »GCON<300) , 
2ttRAD(300) ,VEL(300> »XnPM(1?) ,NKPP(12) ,NCPR(12) , TXK <20. 1 2) . XK ( 20 » 1? ) 

3, TCP (20. 12) .CPX«26,l8) »RHOBX(12).XBM(i2) .FMFBI12) .EMRPil?) ,HXX<1?) 

4 , GAPX ( 12 ) .FTesT (12) .RTEST (12) .TFmDI ( 200) »TXI (POO) ,TX2(200) . 
5TXPTOO.I?) »TUL1 <200J ,TUL2<200».HX(50) .Tw(50l .IPISOI.IRI (501, 

6 TP 2 ( 50) ,TUL(sn) ,IEM<sn) ,TY (200). A<200) »R(200) »C (200) .0(200) , 

?R ( 50 ) , RHO (50 ) ,CP (50) , DXP ( 1 2 ) ,XKR(I0»I?) ,CPB ( 10 , 1? ) . XMOG< 50 ) » 

BYK (50) »AU(10, 12) ,RB( 10. 12), CB (10*12) »DB (10,12), SR (10, 12) , 
9RP1(10,12) ,Rr?(10,12) , H( 1 2 ) ,S (50 ) ,NPM( 12) 

DIMENSION TTUI (50) ,RHOY1 (50) ,RHOY2(50) ,0RH0(50) ,TCPC(20) 

C 

COMMON TKC»XkO» TCPC , CPC , T*V , XX V , TCPV , CPV , XNPM , RHORX , XRM, FMBB , 
1fMFB,NKPb,NCPR,TXK,Xk,TCP,CPX,NPM»GAPX,FTFST»RTEST.TEMDI,TX1, 
PtX?,TX2T,TUL,TUL1,TUL2,IR,IR1.IP2,A,B,C,D,S,R,AP,BB»CP»DB,SR. 

3rP 1 , RR2 , 1 Y . RhoYI , RHOY? , XMD6.RH0 , CP , YK , XKR.CPB , DXR , DT , XL05T , 

4T ABL . TCH AP , TRFC , RHOV , RHOC , FBLOW , FMV , EMC » H300 , NKC , NCPC , NK V , NCPV , 

5nP , NMR , NRPS , nPF , TEST? , TFMP I , TXO , TENV , HEN v , FFNV » OLOSS . TL I M , T I NT 
COMMON II, 12, 13,14,15,15.-0 IN, I NT ,OX,XMT >TL» VL»Bl. .DMP.FRPl .ERR?. 

1 FRR 3 , FRR4 , HV . VPT , CH ARK , CH ARC , ABLK , ABLC • XMDC , H 
C 

CALL PROP 

ynp=np 

S ( I NT ) = ( KHO ( T NT ) *DX*CP (INT )) / ( 2 , 0*OT ) 

R(lNT) = (1.0)/( (DX/2.0)*( (I.O/YKdNT) ) + (1 ,0/YK(INt+1) ) ) ) 

A ( I NT ) SO . 0 

R ( INT ) = ( - ( ( XmOG ( INT ) +XMDC ) *CP ( INT ) +S ( INT ) +R ( JNT )-RHO ( INT ) *CP ( INT) 
1# ( SDOT / ( 2 » 0* ( YNP-1 * 0 ) ) ) ) ) 

C ( I NT ) =XMOG ( I NT+1 ) *CP ( I NT+1 > +P < I NT ) +RHO ( I NT+1 ) *CP ( INT+1 ) *SDOT 
1* ( YNP-1. b/( YNP-1.0 > ) 

D(INT)=<-(GIN+S(INT)*TX?(TNT> ) ) +( XMDG ( INT > -XMDG ( INT+1 ) )*HV 

NPP=NP-1 

JMT=INT+1 

no 10 i=jnt,npp 


*1=1 

S ( I ) = (RHO ( I ) *DX*CP ( I ) ) /DT 

R ( I ) = (l.l))/( (OX/(?,0*YK( I ) ) >+(DX/(2.0*YK(T+l) ) >) 

A < I )=R ( I — 1 ) 

R ( I ) = ( - < XMDG ( T ) *CP (I > +R ( T-l ) +P ( I ) +S < I ) + RHO ( I ) *CP ( I ) *SDOT* ( < YNP-VT 
1-0.5)/(YNP-1.0)> ) > 

C ( I )=XMD6( 1+1 )*CP(I+1 )+P( T)+RHO(I+l)*CP(T+l)*SDOT*( (YNP-Xl-0.5) 

1/ ( YNP-1. 0) ) 

D(I)=(-(S(I) *TX2 ( I ) ) ) + ( XMDG ( I ) -XMDG ( I +1 ) ) *HV 


10 CONTINUE 

R(NP)=(1.0)/( (DXB(l) /(?.0*XKR(l,l|))+<nXB(1) / (2. 0*XKP ( 2 . 1 ) ) )) 

S(NP)=(RHO(NP)*DX*CP(NP)+RHOPX(l )* CPB (1,1) *DYB ( 1 ) ) / ( 2 , 0*DT ) 
A(NP)=R(NP-1) 

R ( NP ) = ( - ( XMDG ( NP ) *CP ( NP) +R ( NP-1 ) +P ( NP ) +S ( NP ) ) ) 

C(NP)=R(NP) 

D(NP)=(-S(NP)+TX?(NP) ) +XMDG(MP)*HV 

no 200 1 = 1 , Nmp 


TF(I-I) 20,20,30 
?0 AP(1,T)=A(NP) 


boooo 

booio 

B0020 

B0030 

BOOUfl 

B0050 

B00B0 

P0070 

P0050 

BOOOO 

P0100 

P011O 

P0120 

P0130 

P0140 

P0150 

P0160 

P0170 

P01B0 

P0190 

P0200 

P0210 

P0220 

P0230 

P0240 

P0250 

P0250 

P02T0 

P02R0 

P0290 

P0300 

P0310 

P0320 

P0330 

P0S40 

P0S50 

P0350 

P0370 

P03AO 

P0390 

P0400 

P0410 

P04?0 

P0430 

P0440 

P0450 

P0460 

P0470 

P0480 

ROUOO 

P0500 

P0510 

P0520 

P0530 

P0540 

P0550 

R0560 
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RR(1,I)=6(NP) P0570 

CP ( 1 » T ) rc ( NP ) P05BO 

nP(l , T)=D(NP> P0590 

fiO TO 65 P0600 

50 i=NPM(I-l) P0610 

IF (FTFST ( T ) ) 45,40,45 P0620 

4 0 BP < 1 , T ) = ( RHORY ( I > *CPR ( 1 , T ) *DXP ( I ) +RHOPX ( T-l ) *CPB ( L > I-t > *0X8 II-IM/ P0630 

1(2. 0*OT ) P0640 

RR1(1»I)=(1.0)/((0XB(T-1 )/(2.0*XkB(L»I- 1 ) ) )+(nXP(I-l)/(2.0*XKP(L-l P065O 

1,1-1)))) P0660 

RP2(l»I)=(1.0)/( (OXB( I)/(?.0*XKP(1 ,1) ) >+{DXR(T)/(2.0*YKR(2»I) ) ) ) P067O 

AP(1,T)=RP1(1,I) P06B0 

RP(1,T)=(-(RR1 ( 1, I >+RR2( 1,1) +58(1,1))) P069O 

CP(1,T)=HP2(1,I) P0700 

0B(1»T)=(-(SR(1,I)*TX?T(1 *1))) P0710 

fiO TO 65 P0720 

45 IF (FTFST ( T ) ) 50,40,55 . P0750 

50 fi=(1.73F-n9)/(1.0/EMRB(T-1 )+1.0/FVFP(T)-1 ,0) P074O 

«0 TO 60 P0750 

55 OSO.O P07AO 

60 BP ( 1 » T ) = ( RHORX ( I ) *CPR ( 1 , T ) *DXP ( I ) ) / ( 2 . 0*OT ) P0770 

RP2(l,I)=(1.0)/( <nXB(T)/(2.0*XKP(l »I) ) > + ( OXB ( T ) / ( 2. 0*YKP < 2 , T ) ) ) ) P07BO 

AP < l , I ) =H ( I-l ) +4 . 0*6* ( TX2T (L , T-l ) **3) P0790 

RP ( 1 , T)=(-(H( T-l >+4,0*6* <TXPT( 1,1 )**3>+RB2(1 , I )+SB( 1*1))) POBOO 

CP(1,T)=RP2(1,I) P0B1O 

nP(l,T)=3.0*G*< (TX2T(|_,I-1 )**4)-(TX2T(l, T)**4) )-RB ( l , I ) *TX2T ( 1 ♦ I ) P0B20 

65 |F=NPM(I)-1 P0B3O 

00 100 J=2»LP POB40 

BP(J,T)=(RHOBY(I)*CPR(J,T)*DXP(I) )/OT POPS" 

RRl(J,I)=(l,0)/( (nXB(T)/(?.0*XKB(J-l»I) ) )+(t)XP(I)/(2.0*XKP(J,T) ) ) ) ROBftO 
RP2( J,I)=(1.0)/< (OXB(I)/ (?,0*XKB( J+lf I) ) )+(DXP(I)/(2.0 *XkR( J,I) ) ) ) P0B70 

AP(J,T)=KB1(J,I) POBgO 

RP< J, T )-(— (R8t (J, T ) +RB2( J, I ) +5B( J, I ) ) ) POB90 

CP(J»T)=RP2U,I> P0900 

DP(J,f >=(-<SR(J,I)*TX2T< J,I> ) ) P0910 

100 CONTINUE P0P20 

IF(I-NMP) 110,250,250 P0930 

110 LNF=NPM(I) P0940 

TF(RTFSTII)) 120,115,120 P0950 

115 BP ( LNF , I ) = ( RHOBX ( I ) *CPB ( |_NF , I ) *OxP ( I ) +RHO«X (1*1) *CPB (1,1+1) *DXB ( T + R0960 

11 ) )/(?.0*DT> P0970 

RPKLNF, I) = (1.0)/( (DXP( I)/(2.0*XkB(LNF- 1 ,T ) ) ) + (DxB ( I ) / ( 2 . 0*XKP (LNF POPfln 

1 , T ) ) ) ) P099O 

RP 2 (LNF,I)S( i , 0 ) / ( (OYP(I+1 >/(2.0*XKB(l »I + l > ) )+(OxB( 1+1 )/ B 1000 

1(2. 0*XKR (2,1+1 ) ) ) ) P1010 

AR(LNF,I)=RBt(LNF,I) P1020 

RP(LNF,I)=(-(RB1(LNF,T)+RP2(LNF,T)+SB(LNF,I) ) ) P1030 

CP(LNF,I)=RB5>(LNF,I) P1040 

(>P (LNF , I ) r ( — ( BB ( LNF » T ) *TX2T (LNF , T ) ) ) P1050 

fiO TO 200 P1060 

120 IF ( RTFST ( I ) ) 125,115,127 P1070 

125 Rr(1.73F-09)/(1.0/EMRP(I)+1.0/EMFB(I+l)-1.0) P10B0 

fiO TO 130 P1090 

127 R— 0.0 R1100 

130 BP(LNF,I ) = (RhOBX( I )*CPB<| NF,I)*DxR(I) >/( 2 . 0 *DT) P1110 

RP1(LNF,I)=(1.0)/( (DxR(I)/(2.0*XKB(LNF-l,T) ))+(DxB(I)/(2.0*XKP(LNF PI 120 
1 , T ) ) ) ) P1130 
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ftP(LNF,!)=RRt (LNFtl) 

RP (LNF » T ):(»(PB1( LNF »T)+H(I)+FB(I. NF*I) +4.0*6* ( TXpT ( LNF » I ) **3 ) ) ) 
CP(LNF f n=H(n+4.n*6*(TX?T(l.T+l >**3) 

OP(LNF»T )=3.fl*G*( <TX?T(1 ,T+1)**4)-(TX2T(I NF,I)**4) >-RB<LNF. I >*TX?T 
l(LNF.T) 

200 CONT1MUF 
250 MN=NpM(NpR) 

T p ( QL0S5 ) 27o,260»27n 

240 <;P < MN » NMB l = ( RHOPX ( NMr ) *CPP ( MN , NMR ) *DXP (NmP )) / ( 2 , n*DT ) 

PP1(MN,NMR)=(1 ,0)/( (0XB<WM6)/(2.*XKP(MN-1 ,NMB) ) )+ (OXR ( NMB) / ( 2. 0*XK 
1r<MN,NMP) ) ) ) 

AP(MNfNMB)=PR1 (MN*NMr) 

RP(MN,NMb)=(-(RP1(MN,NMP)+SB(MN,NMB) ) ) 

CP(MN.NMB)=0.n 

nP(MN,NMfci)S(-(SP(MN»NMB)*TX2T(MN,NMP) ) ) 

GO TO 240 

270 «;P ( MN » NMB ) S ( RWOPX < NMr ) *CPP ( MN »NMR ) *OXR (NmP ) ) / ( 2 , fl*OT ) 

RP1 <MM»NMR)=< 1 ,0>/( (OXB<NMB)/(2,n*XKBCMN-l * NMP ) ) ) ♦ <DXP(N mR) / < 2. 0*X 
lkP(MN'NMb) ) ) ) 

AP(MN,NMB)=RR1 (MN»NMR) 

rP ( MN , NM b ) = ( - ( RB l ( MN , NMP ) +HFNV+ ( 1 » 73E-09 ) *FFNV*4 . 0* ( TX2T ( MN , NMB ) ** 
l.M+FBtMN.NMB) ) ) 
rP(MN»NMri)=0.0 

DP ( MN , NMb ) = < _ < HFNVoTfN V+FFNV* ( 1 , 73E-09 ) * < ( TFNV*t*4 ) +3 . 0* (TX2T ( MN * NM 
1R)**4) )+SP(MN»NMB)*TX2T<MN»NMP) ) ) 

2«0 | rNP+1 

no 30 n 1 =) » Nmp 
K=NPM(I) 

TF(l.c-Q.X) GO TO 282 
TF(GAPX(I-1) .FG.O. ) fiO TO 2B2 
KT=1 

GO TO 2BS 
2«2 K T=2 

245 no 290 J=KT»k 
A(L)=AB( J»I) 

R(L)=«B( J#I) 
r(L)=CB(J»I> 
n(U=nB( j»I) 

TF(DMP) *49»?p9»246 

2«6 *PlTE (6 * 2P7 ) AB ( J* I ) , PB ( J* I ) , rB( J* I ) » OB < I ) . J* I » A <L) »B(L ) #C (l ) *n ( 
11 ) .L 

247 FORMAT(lHO»8HAB(J#I)-,lPF12.5.2X,BHPB(J»ns,lPEl2.5»2X»BHCB(J,I)=» 
llPF12.5»2X»8HnB(J*l)=,lPFl2.5,2X,2HJx»I3.2X»2Hl=,I3/lX#fiHA(L)=tlPF 
2l2.5.2X»5HB(L)=f 1PE12, 5i 2X*5HC (L > =* 1PF12. 5# 2X»5Hn < L )= » 1PF12. 5# 2X.2 
3 hI =»I3) 

249 l =L+1 
2°0 CONTINUE 
300 CONTINUE 
NPFT=L-1 
PFTURN 
PND 


R1140 
P11B0 
BUBO 
PI 170 
R11B0 
R1190 
P1200 
B1210 
P1220 
P1230 
P1240 
B1250 
P1260 
P1270 
B12B0 
P1290 
P1300 
P1310 
P1320 
PI 330 
P1340 
P1350 
P13B0 
P1370 
P13B0 
P1S90 
P1400 
P1410 
P1420 
P1430 
P1440 
B1450 
P1460 
P1470 
P14R0 
R1490 
P 1 5 0 0 
P1510 
P1520 
P1530 
B1540 
P15S0 
PIBfiO 
P1570 
P15B0 
P1S90 
PlfiOO 
PlfilO 
P1620 
P1630 
P1640 
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SltJFTC PRP fOOOft 

this cun * outtme dftfrminfs thf physical properties of thf cooio 

HEAT SHTELD STRUCTURE C0020 

FUBROUTINF PROP C0030 

C C0040 

nlMFNFION TlTt E < 12 ) »HFADNG(12),XIDNT (12,12) »TKC<?0) »XKC(20) . C0050 

lrPC(20) ,TKV<20) ,XKV(?0) ,TCPV(20) ,CPV(20) .TIME (30(1) *QCON(300) , C0060 

2oPAD(300> ,VFl.(300 ) »XnPM( 12),NKPR(12),NCPB(12) ,TXK(20,12) ,XK(20,i?) C0070 

3, TcP(20»12) ,CPX(20, 12) ,RHOBX(12> ,XBM( 12) .FMFB(IP) »EMBP(12) ,HXX( 12 ) COOflO 

4, PAPX(12> ,FTfFT(12> »RTEST(12) .TEMDK200) ,TX1 (200) ,TX?<200). COOBO 

5TX2T(10,12> >TUL1 (200) ,TUl.2(200> ,HX(50) ,Tw(50) , IR(50) , TR1 (50) , C0100 

6lP2<50) , TUL(50) ,IEM(60) , TY < 200 ) » A < 200 ) ,R ( 200 ) ,C ( 200 ) , 0( 200 ) , C0110 

7R ( 50 ) »RHO(50) ,CP(50) ,nXP (12) ,XKB( 10,12) »CPB ( 10 , 1 2 ) » XMD6( 50 ) » C 0120 

RYK(50)»AH(10.12)iRB(10,I2)»CB(10,12)»DB(10*12)»SB(10,12). 00130 

RRP1 (10,12) ,RR2(l0»12) » H ( 12) »S(50) , NPM ( 12 ) C0140 

DIMENSION TTULJ50) »RH0Y1 ( 50 > ,RH0Y2 ( 50 ) ,DRH0 ( 50 ) .TCPC ( 20 ) C0150 

C C0160 

COMMON TKC , XKC » TCPC , CPC , TK V » XK V » TCPV . CPV , XNPM * RHOPX , XPM » FPBP , CO 1 70 

IfPFR , NKPB , NCPP . TXK , Xk , TCP ,CPX ,NPM, GAPX .FTFST , RTE sT . TEMDI , TX1 . C01 BO 

2 t* 2 »TX 2 T,TUL,TUI 1 .TUL 2 , IR,IRl,IR 2 »A»B,C»n,S»R.AP,Re»CP»DR»SB, OOIPO 

3rP1,RP2»TY,RhOY 1 ,RH0Y2,XMnG,RH0,CP»YK,XKR,CPB,DXR.DT ,VLOST, C0200 

4T ABL , TCHAR , TRFC » RHOV/ , RHOC » FRLOW » FMV , EMC » H300 . NKC , NCPC * NK V , NCPV • COR 1 0 

5NP » NMP * NPPS » NPF » TEST 2 » TEMPI , TXO » TENV , HENv , FFNV » Ol.OSS , TLI M » T I NT C0220 

COMMON I 1,12, T3»I4*I5» 16. OIN, TNT. DX,XMT,TL»VL»BL.nMP,FRRt»ERR2, C0230 

1FRR3 ,FRR4 ,HV * VPT • CHARK , CHARC , ABL.K , APLC » XMPC » H C0240 

C C0250 

xTNTrTNT C0260 

HO 170 I=KINT»NP C0270 

10 TF(IR(I)> 12,12,100 C02BO 

12 TUL(1)=AMAX1(TX1(I)»TX2(T)) C02B0 

TF(TUL(I) .LF.TAPL) 60 TO 20 C0300 

TP(I)=1 C0310 

ftO TO 100 C0320 

20 TF(Il-l) 25,25,21 C0330 

21 TF(Il-NKV) 22,22,25 C034O 

22 TF(TX2(I)-TKV(I1)) 35,55.30 C0350 

25 WRITE (6,26) TX2(I) C0360 

26 FORMAT (1H0.P7M THF R ft NGF OF ONE OF THF ARl.ATION PPOPFPTY CURVF FtT C0370 

IF WAS EXCEEDED AT A TFMPFRATURE OF ,1PE12.5) C03B0 

FRR2=1.0 C03S0 

ftO TO 356 C0400 

30 T trll + l C04 1 0 

ft 0 TO 21 C0420 

35 TF(TX?(I)-TKV(I1-1) ) 40,55,50 C0430 

40 T 3 =11—1 C0440 

ftO TO 20 C045O 

60 YK (I )=XKV(I1-1 )♦( (XKV( II >-XKV(Il-l ) ) / ( TK V ( I 1 )-TKv ( I 1-1 ) ) ) C0460 

1*<TX2(I )-TKV( Tl-1 > > C0470 

ftO TO 60 C04B0 

55 YK(I)=XKV(I1) C04S0 

60 IF ( 12-1 ) 25,25,61 C0500 

61 IF ( I2—NCPV) 62,62,25 C0510 

62 TF (TX2 ( I ) —TCPV ( 12 ) ) 70,P5,65 C0520 

65 T2=I2+1 C0530 

ftO TO 61 C0540 

70 IF ( TX2 ( I )— TCPV ( I 2—1 ) ) 75,B5,B0 C0550 

75 T2SI2-1 C0560 
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GO TO 60 

*0 CP ( I ) =CPv ( 12-1 ) ♦ ( ( CPV ( I?) -CPV ( 12-1 > )/( TCPV (12) -TCPV ( 12-1 ) ) ) 

1 * ( TX2 ( I ) -TCPV ( 12— 1 ) ) 

GO TO 90 
05 CP(1)=CPV(12) 

OQ wHO(I)=RHOV 
GO TO 17u 

100 TtlL ( I )=AMAX1 (TUH I) ,TX2< T ) ) 

1F<TUL(I)-TCHAR) 110,110,115 

1 1 0 RHO ( 1 ) =RHOV+ ( RHOV-RHoC ) * ( ( TUL ( I ) -T ABL ) / ( T ABL-TCHAR ) ) 

YK ( I ) =CHARK+ ( ABLK-CH ARK ) * < ( RHO ( I ) -RHOC ) / ( RHOV-RHOC ) ) 

CP ( I ) =CHARC+ ( ABLC-CHARC ) * ( ( RHO ( I ) -RHOC ) / ( RHOV-RHoC ) ) 

GO TO 170 

lib lF(VPT> 116,1)6,117 

116 TTULC T)=TUL(T) 

GO TO 120 

117 TTuLI T)= rX2(I) 

120 IF ( I 3—1 ) 25,25,121 

121 tF ( 13— NKC) 122, 122,25 

122 IF < TH IL ( 1 >— TKC ( 13 ) ) 124,1*5*123 

123 13=13+1 
GO TO 121 

124 TF<TTML<1)-TKC(I3-1)) 125,135,130 

125 13=13-1 
GO TO 12u 

1*0 YK ( I )=XKC( 13-1 )+( ( XKc ( 13 ) -XKC ( 13-1 ) > / < TKC < IS) -TKC ( 13-1 ) ) ) 

1* < TTUL < I ) — TKC ( 13— t ) ) 

GO TO 140 
135 yK(I)=XK0(I3) 

140 IF ( I 4—1 ) 25,25, )4l 

141 IF ( T 4— NCPC ) 142,142,25 

142 lF(TTltL(l)-TCPC(I4) ) 150,165,145 
14b 15=14+1 

GO TO 141 

150 )F (TTUL ( 1 )— TCPC ( 14—1 ) ) 155,165,160 
155 14=14-1 

GO TO 140 

160 rP{I>=CPC<I5-i >+( <CPC( 15 )-CPC( 14-1))/ (TCPC <I4)-TCPC< 14-1 ) ) ) 

1 * ( TTUL ( I ) -TCPC ( 14-1 ) ) 

GO TO 166 

165 C.P(I)=CPC(I4) 

166 WHO ( I ) =PHOC 
170 CONTINUE 

nFTFRMlNATIOiM OF PROPFR BACK-UP 5HIFLP MATERIAL PROPERTY 
no 300 1=1, NMR 

nXB ( I ) =XBM (!)/(( XNPM ( T ) — 1 .0 )*12.0) 

LKP=NKPB(I) 

I CP=NCPP(I) 

nn=npm< i ) 
no 260 J= 1 ,Ni\| 

200 iF(Ib-l) 203,203,201 

201 1F( Ib-LKP) 202»2n2»2n3 

202 1F(TX2T(J,I)-TXK( 15,1) ) 206,220,205 

203 ivPITE ( 6,204 ) T,TX2T(J,I) 

204 P ORMAT ( 1 HO , 32H THF RANGE OF ONE OF THF NllMBFR ,I2»71H BACKUP 5TPUC 


C0570 

C05B0 

C0591) 

C0600 

C0610 

C0620 

C0630 

C0640 

C0650 

C0660 

C0670 

C06B0 

C0690 

C0700 

C0710 

C0720 

C0730 

C0740 

C0750 

C0760 

C0770 

007*0 

C0790 

COBOO 

CORK) 

00*20 

00*30 

COR40 

COB50 

C0B6O 

C0R7O 

00**0 

COB90 

COROO 

C0R10 

C092O 

00930 

C0940 

00950 

00960 

00970 

CORRO 

C0R9O 

oiono 
cim n 

C1020 
C1030 
C1040 
C1050 
C1060 
C1070 
C10R0 
C1090 
C1100 
C1110 
Cl 120 
C 1 130 
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1TURF PROPERTY CURVE FITS WAS FXCFFDFD AT A TFMPFpATURF OF .1PF12.5 01160 

2) 01150 

FRR2=1.0 01160 

fiO TO 35b 0117(1 

2(15 Tbslb+l 01180 

fiO TO 201 01190 

206 TFUX2T<J,I)-TXK(T5-1,I) > 210,220.215 01200 

210 T5=Ib-l C1210 

fiO TO 200 01220 

215 XKR(J,I)=XK<lF-l,I> + ( (XK(T5»I)-XKU5-1,I) )/<TXK<T5.I)-TXMl6-l,T> ) 01230 

1 >* <TX2T<d, I >-TXK< 15-1,1) ) Cl2ttfl 

GO TO 230 01250 

220 xKBCJ,I)=XK<I6,I) 01260 

230 TF(16-1) 203,203,231 01270 

231 TF(Ib-LCP) 232,232,203 01260 

232 TF(TX2T(d,X)-TCP(Ib,n) 234,245,233 01290 

233 Tb=lb+1 01300 

r,0 TO 231 01310 

234 TF(TX2T(J,I)-TCP(Ib-1 ,11) 235.246.240 01320 

235 T6=I6-1 01330 

fiO TO 23U 01340 

260 CPB(J,I)=OPX( T6-1 ,!)♦( <PPX<I6»I)-CPX(T6-1 »I) ) / ( TOP ( 16 » I )-TCP ( T6-1 . 01350 

IT) ) >*<TX2T<J,T>-TCP<T6-1,T>> 01360 

fiO TO 260 01370 

245 <*PB(J,1)=CPX(T6»I> 01360 

260 CONTINUF 01390 

T5=2 01400 

T6=2 01610 

300 CONTJMUF 01620 

310 TF(DMP> 365,365.320 0 1630 

320 wPlTE<6»330) 01640 

330 FORMAT ( /IX , 32M PROPFRTIFS OF ABLATION MATFRTAL/) 01650 

wPITE(6,335) 01660 

335 FORMAT (/5X,5HYK( I) »9x,5H0P(I) .9X.6HRH0 C I ) / ) 01470 

WRITE (6,340 ) (YK(I),CP(I).RH0(I),I=1,NP) 01460 

360 F0RMAT(2X»1Pf12»5.2X,1PF12»5»2X,1 PEI 2. 5) 01490 

WPITE(6,365) 01500 

365 F0RMAT(//1X,32H PROPFRTIFS OF BACK-UP STRUCTURE/) 01610 

wPITE (6. 367 ) 01520 

367 FORMAT ( /5X , 6hyKB ( J, I ) , 7X , AHCPR ( J, I ) .7Y .RhRHOBX ( I ) ,7X , THFmFB ( T) » 6x » 01530 

17HFMBR< I) ,9X,6HDXR< I) /) 01560 

00 350 1=1, Nmb 01650 

Kl =NPM ( I ) 01560 

no 349 J= 1 ,KL 01570 

WRITE (6, 368) YKP ( J, I ) ,CPR < J, I ) .RHOBX ( T ) ,FMFB ( T ) ,FMBR (T ) ,DXB < I ) 01660 

368 F0RMAT(3X,1PF12.5,3X,1PE12.5,3X,1PE12.5»3X,1PF12,R,3X,1PF12.5,3X.1 01590 

1PF12.6) 01600 

369 CONTINUE 01610 

360 CONTINUF 01620 

355 RETURN 01630 

FMP 01640 
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SIBPTC ABL 

THIS SUBROUTINE DETERMINES THF MASS FI OW PATF FROM THF 
ABLATING NOnES 
SUBROUTINE ABI ATE 

DIMENSION TITLE ( 1? ) ,HFADNG ( 12 ) »X tDNT (12,12),TKC(20),XKC(20) , 
lrPC(20) ,TKV(20> ,XKV(20> ,TCPV(20).CPV(20) , TIME <30O)»QCON( 300) . 
2(lPAD<300) ,VEL(300> »XNPM(12) »NKPB(12> »NCPB(12> ,TXk(20, 12) ,XK(20, 12) 

3, TCP{20»12> ,CPX(20»l2> ,RH0BX(12> ,XBM(12> .FMFBC1P) ,EMBB(12) ,HXX(l2) 

4 , GAPX ( 1? ) »FTfsT ( 12 ) »RTE5T ( 1? ) ,TEmOI ( 200 ) , TX1 ( 200) * TX2 I 200 1 . 
5TX?T(10»12> ,TuLl(200) ,TU| 2(200) ,HX(50) ,Tw( 50) , IR( 50 ) , TR1 ( 50) , 

6tP 2 ( 50 ) , TUL ( 50 ) » IEM(sf)) ,TY(200) , A(200) ,B(200> ,C(p00) ,O(2n0) , 

7r(50> , RHO (50) , CPI 50) ,0XP(12> ,XKB< 10,12) ,CPB ( 10 , 12) • XM06I 50) * 
BYK(50> .AtM 10,1 2) ,RB( 10.12 ) »CR( 10. 12).0B( 10.12 ) ,SR( 10,12) , 

9RP1 (10,12) .RR2I10.12) ,H(12) *5(50) »NPM( 1?) 

DImFNSION TTULI50) *RhOY1 I 50 ) ,RHOyP ( 50 ) »0RH0 ( 50 ) ,TCPC < 20 > 

COMMON TKC * X«C * TePC . CPC , TK V * XK V , TCPV * CPV , XNPM , RHoRX , XBM * F M BR , 
1FMFB,NKPB.NCPB,TXK,Xk,TCP*CPX,NPM»6APX.FTFST.BTEST*TEMDI.TX1, 
2TV2.TX2T ,TUL, TULl , TUL2* TR»IR1»IR2»A»B,C,0,S»R, AB.BB *CB ,DB* SB, 
3 rB 1,RR2,TY, RH0Y1 . RHO Y2 » XMOG , RHO , CP » YK , XK B , CPB , DXR , DT , XLOST * 
4TABL,TCHAR»TRFC»RHOV,PHOC.FRLOW,FMV»EMC,H300»NKC,NCPe.NKV,NCPV, 

5NP , NMB , NPBS * NPF . TFST2 , TFMPI * TXO . TENV, HENv , FFNV , OLOSS , TLIM, TINT 
COMMON 1 1,12, 13.14,15, IS. OIN, JNT . OX, XMT »TL. VL .BL.OMP.FRRl .ERR2. 

IF RR3 , FRR4 » HV , VPT , CH ARK , CM ARC . ABLk * ABLC , XMO« ,JH 

XMT=0.0 
l.TNTsTNT 
kI=NP 

T F ( DMP ) rt » 8 , 3 
3 WRITE (6.5) 

5 FORMAT (//IX, goHMASS FLOW FROM ABLATING NODES//) 

8 no POn KKT=LlNT,NP 
TF(IRKKI)) 11.11,12 

11 TF(TX1 ( K I ) .LF.TABL) 60 TO 9 

12 TULl IKI )=AMAX1 (TULI(KI).TYI(KI)) 

TPllKDsi 
60 TO 20 

9 TFITX1 (K 1 )-TABL) 10,10,20 
10 RH0Y1 (KI )=RHOV 

GO TO 50 

20 TF(TULKKI)-TCHAR) 40, 30, SO 
30 RHOYKKI )=RH0C 
GO TO 50 

40 RHOY1 (KT )=RHOV+(RHOV-RHOC)*( (TULl (KI)-TA rI )/(TAB| -TCHAR) ) 

50 TF(IR?(Ki)> 52,52,54 
52 IF(TX2(KI),LE.TABL) GO TO 56 
54 TUL2(KI)=AMAX1 (T|JL2 (kI) ,TX2(KI) ) 

IP2(KT)=1 
GO TO 70 

56 TF(TX2(KI)-TABL) 60,60,70 
60 RHOY2IKI )=RH0V 
GO TO 95 

70 TF(TUL2(KI)-TCHAR) 9n,80,«0 
BO RHOY2IKI )=RHOC 
GO TO 95 

<30 RH0Y2 ( K I ) =RHOV+ ( RHOV-RHOC ) * ( ( TUL2 ( K I ) -T ARl ) / ( T ABl.-TCHAR ) ) 


nooon 

nooin 

D0020 

D0030 

00040 

noosn 

D0060 

00070 

ooobo 

D0090 

noinn 
noun 
00120 
no iso 
noiun 
D0150 

noi 6n 

D0170 
00180 
noi 90 
D0200 

no 2 i n 

D0220 
n0230 
D0240 
00250 
D0260 
D0270 
D02B0 
00290 
D0300 
no3i n 
D0320 
no 330 
D0340 
D0350 
D0360 
D0370 

nosan 

00390 

nouon 

00410 

00420 

00430 

00440 

00450 

00460 

00470 

00480 

00490 

00500 

00510 

00520 

00530 

00540 

00550 

00560 
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95 nRHO(KI) = C (RhOY1(KI>-RH0Y!><KI> >/OT)*D* 00570 

TF(KI-NP) 97,06.96 00580 

96 nPHO(KI)=ORHn(KI>/2.0 00590 

60 TO 9R 00600 

97 TR(KI-INT) 96,96,98 00610 

98 IP (DRHO(KI ) ) 110,120,120 O06?O 

110 ORHO(KI)=0.0 00630 

120 XMT=XMT+ORHO(Kl> 00640 

XMDG(KI>=XMT 00650 

TF(DMP) 190 » 190 ,150 00660 

150 wPITE (6» 160 ) yM06(KI),DRMO(KI)*RH0Y2(Kn,PH0Y1(KT) 00670 

160 FORMAT (IX, 5HXMDG= , lPF 1 2 . 5 , 2X , 5HDRH0= , 1 PF 1 2 . 5 » 2X » 6HRH0Y2= , 1 PF 1 2 . 5 . 2 00680 

1X,6HRM0Y1=,1PF12.5) 00690 

190 Kl=KI-l 00700 

200 CONTINUF 00710 

ftFTURN 00720 

phiD 00730 
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IBPTC OX ID FOOOO 

room 

THIS SUBROUTINE CALCULATES THF HFATING RATE DUE TO COMBUSTION FOORO 

IT IS ASSUMED THAT OXYGEN AND CARBON REACT TO FORM CO ONLY. F0030 

F0040 

SUBROUTINE OXTDAT(XMDO,OOXID) F0050 

C FOOfiO 

ooxiD=XMuo*4ono,n/3Fnn.n roo 7 n 

OOXIDrO.O FOOan 

RFTURN F0090 

END F0100 
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4IBFTC RwUFT 

THIS BUPROUTIME DFTFrMINfF THF FoPWARP TT«E FTFP TEMPFRATURFR 
RY SOLVING ThF TRI-DiAGONAL MATPTX 
BURRQI )T I NF SWUFT ( A , R , r r P , T » N , PMP ) 

dimension A(?no> *R(?nn) ,c<200> .D(?on> ,T(?oo> ,rp(?no) »nP(?no> 
CPU)=C(1)/R(1 ) 
nP«l)=o<i)/R(i) 
no 100 T=?*N 

rP( I >=C< 1)/(R( I )-A( I )*CP< T-l ) > 
nP(i) = (P(T)-A(l)»nP(r-i> )/<R<n-A(l)*cP<T-i) > 

100 rONTIMUF 
T(N)=PP(N) 
mM1=N-1 

no 200 j=i f nmi 
1=N-J 

T I I >=PP(i )~CP ( I )*T ( 1 + 1 ) 

200 ('ONTINUF 

TF(PMP) 300»300»250 
250 wPITE ( 6 » 2f>0 ) 

2F0 FORMAT(//l Xm^HCOFFFTriFNTS CALCULATED By SUBROUTTNF SWUFT//) 
wR ITt ( 6 » 270 ) 

270 FORMAT ( AX , 5HCP ( I > » 1 0 X . 5HPP ( 1 1 , 1 0 X # 4HT ( I ) /) 
wPITE ( 6 » 275 ) (CPU) ,nP (I) .T <I> . 1=1 iN) 

275 F0RMAT(2X»1PE12.5*2X, 1PF1P.5#2X» 1PE12.5) 

300 PFTURM 
END 


F0000 

F0010 

F0020 

F003O 
F00UO 
F00RO 
F00AO 
F0070 
FOOftO 
FOOPO 
F010O 
FOUR 
F012O 
F01 30 
F0140 
F0150 
FOlfiO 
F0170 
F01 AO 
F01PO 
F020P 
F0210 
F0220 
F0230 
F0240 
F0250 
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ibFTc pec oonno 

oooio 

this subroutine determines the front facf location a no char mass 00020 

rfmoval hate 00030 

00040 

SUBROUTINE RFCESS(XMnC»X| OST *TRFC* DT ,RHOC *TS» SR »TX2»NREC , NRS *FRRs * 00050 

lSXO»SDOT,nMP) 00060 

C 00070 

DIMENSION TS(kO) ,SR(50> 00080 

TF(TX?-TKEC) 10»?0»20 00090 

10 XMDCSO.O OOIOO 

XL OST =0 • II 80110 

SDOTsO » 0 00120 

RO TO 60 00130 

20 IF (NRS-1 )2b»25»2l OOlUO 

21 TF(NRS-NkEC) P2 fPPtP* 00150 

22 TF(TX2— TS»(NRS1 ) 32»4(l»30 00160 

25 wPITE (6*26) TX2 00170 

26 FORMAT ( 1H0 * 7SH THE RANGF OF THE SURFACE RFCFSSION TARI E WAS FXCFfO 00180 

1FD AT A TEMPERATURE OF , 1PE12.5) 00190 

FPR5=1.0 00200 

00 TO 60 00210 

30 NPS=NRS+1 00220 

00 TO 21 00230 

32 TE(TX?-TS(NRs-1) ) 34,40*36 80240 

34 NPS=NRS-3 00250 

RO TO 20 00260 

36 5X=SR< NRS-1 )♦< < SR < NRsi -SR ( NRS-1 >>/< TS ( NRS) -TS < NRs-1 > >1 00270 

1 *( TX2—TS ( NRS-1 ) ) 80280 

rO TO 50 00290 

40 SXrSR(NRS) 00300 

SO xl OST=300.0*S**DT 00310 

xMDC=(XLOST*RHOC)/DT 00320 

FD0T=SX*300,n 80330 

TF(DMP) 60»60,52 80340 

62 wR ITt (6 » 54 > SV,XL0ST,VMDC 00350 

54 FORMAT (1 HO » 3HSX= , 1 PE 1 2 . 5 , 3X , 6HXL0ST= , 1PE 1 2 . 5 , 3X » 5HXM0CS * 1 PEI 2 . 5) 00360 

60 RFTURN 00370 

FND 00380 
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IBFTC TEMP H0000 

THIS SUBROUTINE DETERMINES THF INITIAL TF-MPFR ATllRF DISTRIBUTION HOfllO 

IN THF HEAT SHIELD STRUCT! tRF H0020 

SUBROUT INF TEMPD • H0030 

C H0040 

DTMFNSION TIT! E(1?) »HF AOnP ( 1 2 ) » X IDNT (12»1?)»TKC(P0) »XXC(?0) » H0050 

ir.pc(2n) * iKV(2n) >xkv(po) ,tcpv(?o) ,cpv(po) .timfooo) .qcon(300) . hoobo 

2oPAD(300> ,VEL<300) .XNPM(IP) *NXPB(12)»NCPR(1?> , TXX ( 20 . 1 2 ) , *K ( 20 , tp 1 H0070 

3tTCP(P0»i?) »CPX <?0rlP)» PHOBX (12)»XBM(12)» FMFB ( 12 ) .EMRP( Jp) ,HXX< 1.9) H00A0 

4 , BAPX ( 12 ) »FTEST ( 1 ? ) » RTEST ( 12 ) tTEMDI (200) .TX1 (POO) »T¥2<200) , H0090 

5TV2T(10»1?) fTULl (200) ,TUt 9 ( 200 ) , HX ( SO ) * Tw { 50 ) , IR ( SO ) , TR!< 50 ) , M0I00 

6lR2 (50 ) » TUL ( 50 ) » IEM( SO) » TY (200 ) > A ( 200) »R(200) »C(200) ,n(?00) , HOI 10 

7r(50) »RHO(50) ,CP(50) ,DXP(12) ,XKB( 10»12) »CPB(10,12) »XMO6(R0) , H0120 

RYK(50) »AW(10.12) »RB(10,12) ,CR(10,12)»OB(10*12) »SR(10.t2). HOI 30 

9RPJ ( 1 0 » 12) »RR2( 10 » 12) # HI 1 9) *S (50 ) » NPM( 12 ) HOI 40 

DIMENSION TTllt ( SO ) >RH0Y1 ( SO ) .PHOY2 ( 50 ) *DRHO ( 50 ) »TCPC ( 20 ) H01SO 

C H0160 

COMMON TkO»XkC» TcRC>CPC»T4V»XXV»TCPV»CPV,XNPM»RHnPX*XPM»FMBB, HOI 70 

1FMRR,nkPh»NCPR»TXK»Xk,TCP»CPX,NPM»GAPX.FTFST.PTFsT»TEM 0I,TX1 , HOlflO 

2tX?,TY2T»TUL,TUL1»TI)L2* IR»IR1 ,IR 2 . A»B*Cfn»S.R»AP,RB#CR#08#SB# HOI 90 

3rR1»RP2»1 Y,RH0Y1 ,RH0Y9»XMnG.RH0»CPfYK*XKR l CPB.DXB,0T,XL0ST# H0200 

4TARL »TCHAR,TRrC»RHOV.RHOC#FBLOW»FMV<EMC»H300.NKC,NCPC»NKv,NCPV/f H0910 

BNP, NMB # NPRS » NPF »TFST?,TFMPI»TXO»TFNV»HENv#FFNV# GLOSS »TLIM» TINT HOPPO 

COMMON II »I2,I3» I4»IS» If>*OlN» TNT #DX* XMT »TL • VL * BL . OMP , FRRl #ERR2» HOP30 

1 F RR3 , FRR4 » HV , VPT , CH Ark . CH ARC » ABLK » ARLC » XMDC » H H0240 

C H02S0 

*=0.0 H0260 

TF(TESTP) 30n,10n.20n H0270 

100 DO 150 L=1»NPF H0280 

TXl(L)=TtMPI H0290 

TXp(L)=TtMPI H0300 

TUL1 (L)=TX1 (L) H0310 

TUL?(L)=1X2(L ) H03P0 

TFMDI (L)=TEMPI H0330 

ISO CONTINUF H0340 

DO IfeO 1=1 »Nmb H03S0 

V |N=NPM ( I ) H0360 

DO 15S M=1 , JN H0370 

TX2T(M»I )=TEMPI H03RO 

1S5 CONTINUE H0390 

ISO CONTINUE H0400 

00 TO 32 U H0410 

200 DO 220 L=lfNP H04P0 

TFMDI (L)=TXO+( (TFNV-TXO)/TL)**X H0430 

TX1(L)=TEMDI(I ) U0440 

TX2(L)=TX1 (L) H04SO 

TIIL1(L)=1X1(L) HOAfiO 

TUL2(L)=TX1(L) H0470 

X=X+DX H04R0 

2?0 CONTINUE H0490 

I =NP+1 H0500 

no 270 I=1»NMR H0S1O 

KJ=NPM(I) H0520 

DO 250 J=1 iKj H0530 

TFMDI (L)=TX0+( (TENV-TXO)/TL)**X H0540 

TV1 (L)=TEMDI (I ) H0550 

T X p ( L ) =T E MD I ( L ) HOSfin 
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TV2T( J»I )=TEMnl(L) H0S70 

X=X+DYB(1) M05AO 

L=L+1 HOMO 

2*0 CONTIMUF HOMO 

X=X+(RAPX(I)/t2.0) M0610 

270 CONTINUF H06?0 

BO TO 320 H0630 

AM ARRRTTARY TEMPFRATURF DISTRIBUTION CAN BF PEAn IN FROM INPUT M06U0 

DATA IF TEST? IS A NF6ATTVF NUMRFR H0650 

300 uiPITE ( 6» 31 0 ) HOAAO 

310 FORMAT < 1 HO »7«*H THF VALUF OF TFST? WAS NFgATTVF# SUBROUTlNF TEMPO S HOA70 

IHOULD NOT HAVF BFFN CALLFn.) HOMO 

FRR1S1.0 HOMO 

3?0 rFTURN H0700 

FMD H0710 
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ilBFTC P0N2 

THIS SUBROUTINE DFTFRMINfS THE TfMPERATUrF OF POINTS A FIXED 
DISTANCE FROM A RFFFRFNCF PLANE FROM THE TEMPFR ATURES CALCULATED 
IN A VARYING THICKNESS 

Si 'BROI ITINF DnN 2 < XLOS T » XApPAY . T ARRAY , NA , XnODF , TFMp , XNonEV. XK , XI STv » 
1DX) 

C 

DIMENSION XARPAYC50 ) .T APR AY (50) »XNODE(50> # TEMP (50) »XNODFV(50 ) 

C 

K=0 

DXT=0.0 
do ton 1 = 1 »na 
rFtXLSTV.LE.DXT) go to 1RO 
K=K + 1 

ino dxt=dxt+ux 

ISO kK=NA-K 

XK=K 

xNODEV ( 1 ) rXLSTV 
TFMP( 1 )=TARRaV( 1 ) 

DO ROD 1=1 rKK 
XNODE 1 1 ) =XK*DY-XLSTV 

CALL DISCT3(XNODE(I) ,XARRAY»TARPAY.NA*TFMP(T + 1 ) ) 

XNODEV (1+1 >=XK*DX 
200 XK=XK+1.0 
WFTURN 
FND 


looon 

loom 

T00R0 

T0030 

10040 

lonsn 

10060 

10070 

TOOflO 

T0090 

T 0100 

T0110 

I0I20 

TD130 

10140 

T015O 

T0160 

10170 

TOlftO 

T0190 

TOROO 

T0R1O 

T0R20 

T0230 

T0R40 

TORSO 

T0260 
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SXBFTC 

UINTRP 

jonno 


6UBROUTINF UIMTRP(X»XTBL ,Y»YTPL»N»J> 

jooin 


niMFNSION XTRl ( 50 ) » YTPL ( 60 ) 

J0020 


I=J 

joosn 


TF(I.6T.N.0R.T.IT.2> t=? 

J0040 

10 

TF ( XTPL ( I— 1) .1 E»X. ANn, X, l F #XTPL ( 1 1 ) 60 TO 40 

JOOPO 


TF<X.6T,XTBL(T) ) 60 TO SO 

jOOfiO 

20 

T = I-1 

J0070 


TF(I.6E.2> 60 TO 10 

joofln 


T -? 

jooqn 


60 TO 40 

joioo 

30 

T-I + l 

jotio 


TF(I.LE.N) 60 TO 10 

JO 120 


t=N 

joiso 

40 

PRACT=(X-XTPL(I-1) )/(XTPL«l)-XTPL(I-l) ) 

joiun 


y=ytbki-i h-( ytpl< d-ytpi <i«i > >*fract 

J0150 


RFTUKN 

JO 160 


FK’D 

J0 1 70 
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SIBFTC 1S0T KOOOO 

m IBROUTINF ISOTHMC DEPTH, TFMP,PONn»N> ANSI 40010 

DIMENSION DFPTH(l) »TFMP(1 1 400P0 

ANS=-1 . 40030 

K=N-1 40040 

no ion t=i»k Koosn 

TF(TEMP(l)-BOND)?,l,3 KOOfin 

1 ANSSDFPTH(I) 40070 

oo to loo KOO«n 

2 lF(TEMP(l+l)-POND)10n,lnn,4 40040 

4 AMS=DFPTH ( I + 1 ) - < TFMP ( I + 1 ) -BOND ) * ( DEPTH ( I + 1 ) -DEPTH < I >> / < TpMP < I + 1 1 - 4 0100 

1TFMP<T)) 40110 

00 TO 100 K0120 

3 TF(TE«P(1+1)-B0ND)5, 100,100 40130 

5 AN*r ( TEMP ( I ) -ROND ) * ( DFPTH < I 1 ) -DfPTH (!))/< TFMP( I ) -TEMPI 1 + 1)1 +DEPTH K 0 1 40 

KI) K01S0 

100 CONTINUF 401 AO 

TF(BOND.eO.TFMP(N> > ANRrDFPTH(N) 40170 

RFTURN 401BO 

PND 40140 
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SIBFTC FAVE L 0000 

FUBROUTINF SAVE(FAVFl ,SA vF2*SAVE3»USE,NX1 ,VALUE*OT*TFTNA| ,TIMF» I 0010 

1THING) 10020 

DTMFNFION SAVF1 « 1 ) »SAVE2 ( 1 ) »SAVE3( 1 1 l.OOSO 

HSE=0.0 l OOttO 

FAVE1(NX1)=VAIUE I.OOBO 

nX2=NX1-1 L0060 

TF(NX?.FQ.0)NX2=3 1.0070 

FA VE2 ( NX2 I'SVAl UE 1.0080 

MX3=NX2-1 l 0090 

TF(NX3.EG.0)NV3=3 10100 

FAVF3(NX3>=VA| UF 10110 

TF( (TTME.IT. <?.*DT) 1 .OR. (TIMF.GF. (TFINAL-3.*DT1 1 )60 TO 4 10120 

GO TO (1 »2»3) ,NX1 10130 

1 TF( ( (»BS(FAVF?(1 >-SAvF 2 (?) 1 1 .I.E..001) .OR. (ABS(SAvE2 (21-SAVE? (3) 1 1.01 40 

1.1 E..001) 1G0 TO 5 101 BO 

TF( < (FAVE?(1 1 ,LT.SAVf?(?> l.ANO. (FAVF2(2) .PT.SAVF?(31 1 1 .OR. ( (FAVF?< I 0180 

111. T,SAvF 2(21 ) . AND, (SAVF?(21 ,LT. SAVE? ( 311)1 USEssAVE? ( 2 1 10170 

b TH1NG=SAVF2(21 LOlftO 

GO TO 4 I 0190 

2 TF< ( (ABS(FAVFX(l)-SAvF3(?) ) ) .1 E. . 001 1 .OR . ( APS ( SAvE3( 2 1-SAVE3 < 3) 1 L0200 

1.1 F..001 1 )G0 TO 6 L0210 

TF ( ( (FAVE3<1> .LT.5AVF3(P) l.ANO. (FAVE3(2> .6T.SAVF.K3) 1 1 .OR . ( ( SAVF3 ( (.0220 

11 1 .PT .SAVF3 ( P 1 ) . AND, ( FAVF3 (21 ,LT.SAVE3<3> 1 1 1USESfAVF3(2> I 0230 

6 THING=SAVF3(?1 L 0240 

GO TO 4 1.0250 

3 TF( ( ( AB5(5AVF1 (ll-SAVFH?) 1) .1, E. . 001 » .OR. ( APStSAvFl (21-SAVE1 <3> ) L0260 

1 .1 E. .001 ) 1 GO TO 6 L0270 

1 F ( ( (FAVE1 (11 .LT.FAVF1 (21). AND. (FAVEK2) .PT.SAVFl (31)1 .OR. ( (SAVE 1 ( L02R0 

111 .PT.SAVFl (?) ) .AND. (FAVF1 (21 .LT.SAVEl (3) 1 1 lURESFAVFl (2) (.0290 

7 TH1NG=SAVF1(2) I 0300 

4 N;XirNXl+l ( 0310 

TF(NXt.F0.4)NVlri 10320 

RFTURN L0330 

END L0340 
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SlbFTC PISCT3 

FtlBROttTTNF DTFCT3 ( xa , tab* »TARY»NY» ANS) 
niMFNFION T ABY ( 1 ) # T ARY ( ) ) 

CALL niSSFR(X*»TARX*l »NY,J>.NN) 

NNN=3 

CALL LAGhAN(XA,TAPX(nN) *TABY(NN) ,NNN. AN?) 

RFTURN 

FMD 


«ooon 

mooio 

MOO?n 

M003rt 

M004n 

MOORfl 

M0060 

M007n 
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SIBFTC MORE NOOOO 

DIMENSION TITLE < 12) » X( 2000) » Yl (2000) »Y2<2000> ,Y3(2000) *Y4<2000) N0010 

RFWlND ll N0020 

HFADCll) (TITLE(I)#I=1»12> NOOSO 

RFAD ( 1 1 ) X (1 ) # Yl ( 1 ) » Y2 ( 1 > , Y3 ( 1 ) » Yu ( 1 > NOOttO 

YS(l)rY3(ll*t?.+YH1l NOOFO 

Y4(1)=Y4(1)*1?.+Y1(1) NOOftO 

T=2 N0070 

■30 RFAD(11)X(I)*Yl(I)»Y?(I> ,Y3(I>»Ytt(I) NOOftfl 

TF(X(T>-**001.)10»20r?0 NOOQO 

10 y3(I)=Y3(I)*1?.+Y1(I> M0100 

YU(I)=Y4(T)*1?.+Y1U) WOliO 

7=1+1 M01 20 

00 TO 30 MO 130 

20 MPL0T=I-1 NOlan 

YM1=Y1(1) NOlfifl 

YMp=Y?(l) NOlftO 

YM3=Y3(1) MO 170 

YM4=Y4(1) NOlflO 

no 40 K = 2 . NPLOT N01RO 

TF (Yl ( K > .6T.YM1) YMi = Y1(K) M0200 

TF (Y?<K ) .GT.YM2) YM? = Y?(X) M0210 

TF (YMK ) .6T.YM3) YM.3 = YMK) M0220 

TF (Y4(K) .GT.YM4) YMu = Y4(K) M0230 

40 CONTINUE M0240 

1000 F0RMATUH1 » (l?Af ) ) N02S0 

CALL ACCEND (X»Y1»Y2»Y3» Y4 * NPLOT ) N0260 

XMftX=Y (NPLOT) M0270 

CALL APLOT (X»Y1, XMftY, YMI .TITLF, NPLOT) M02A0 

CALL PPLOT <X,Y2,XMAy, YM p, TITLE) N02RO 

CALL CPLOT ( X » Y3» Y4 » XMAX * YM3» YM4 » TITLF» Yl ) M0300 

WRITE (6» 1000) (TITLE(T) f 1=1 »1?) N0310 

WRITE (6 » 1001) (X(I) #Yl (I) ,Y2(I) *Y3(I) *Y4( T) * 1=1 ♦ NPLOT) N0320 

1001 FORMAT.! 5E20 , A) N0330 

WRITE (6» 1002) XMAX»YM t, YM? * Y M3 » YMu » NPLOT N03U0 

1002 FORMAT ( ///6H XMAX=Fl 0.4*FH YM1=F10.4,RH YM2=F10.4»SH YM3 sF10.4*5H N03S0 

lYM4=Fl0.4»2Xf>HNPLOT=I4) N03A0 

RFAD (11) (TITLE (I),T = 1 » 12 ) N0370 

rFAD( 11)X(1) »Yl( 1) »Y2( 1) ,Y3(1 ) »Ya(1) M0380 

T=2 N03RO 

TF(X(1)-S001.)30»50»50 N0400 

FO wPITE(6»1003) ( TITLE ( T ) » 1 = 1 » 12 ) N0410 

1003 FORMAT (////12A6) N0420 

rFTURM N0430 

FMD N0440 
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SIBFTC ACCEN POOftft 

fubroutinf accend<x»y»a»b»c»ni ponin 

nlMENFION X(1WV(1)»A(1),PC1).C(1) POO?ft 

K=1 POOSft 

101 <;MALL=X<K) POOuO 

no ion i=k*n POOsn 

nuMYsxd) pooftn 

FMALLsAMjNl (SMALL »OU mY) P0070 

tf(Small.fq.x(I))Inofysi pooan 

ino CONTINUE P0090 

X(INDFX)=X(K) POIOO 

x<k>=cmali. poiin 

faveey <k ) potsn 

Y ( K )SSY ( INDEX ) P013n 

yundfx)=«;ave pomn 

«;AVFA=A(k) Pflisn 

A(K)£A(lNOEX) poiao 

A(INDFX)sFAVFA P017O 

*;AVFB=B(K) POlAft 

R ( K ) SR ( INDEX ) P0190 

R(INDFX)=FAVER P0200 

qavfcscik) posin 

c(K)=r(tNnEX) po??n 

CCINDFX)sFAVEC P02SO 

KSK+1 PO?ftft 

tF(K.fq,n)Retmrn PO?*n 

GO TO 101 P02AO 

FMD P0270 
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SIBPTC APLOT 

SUBROUTINE API OT ( X , Y , XL TM, YLTM»TITLE» IPl OT ) 

DIMENSION X(300> .YTlTLE(tO) ,XTITLF(10) 
nTMFNSION TITIE(12>»Y(300>»ALONGY(7) 

COMMON /ARC / ALLOW (7) » Al ONGX (7 ) »NPLOT » ?FRO» XMAX, IFIX 

nATA <XTITLE(T)*I=1»10)/3RH TIME (SEC.) / 

nATA (YTiTLFn) ,I = 1»10)/38H SURFACE RFCF SSI On (IN.) / 

?FRO=0.0 

Al L0w(l>=50. 

AlLOw(2>=!00. 

Al L0w<3>=250. 

Al LOW (4) =500. 

AILOw(5)=100o. 

At LOW (6)=2500. 

Al LOW ( 7 ) =5000 . 

NPLOT=IPLOT 
DO 10 1=1 >7 
T 1 = 1 

TFIXLTM- ALLOW ( I ) ) 2n,20,10 
10 OONTINUF 

30 wPITE (6,1000) XLIM,YI.IM 

1000 FORMAT ( ///77H APLOT CANNOT FF DONF BECAIJSF FITHFr XLTM FXCEFOFO 50 
Inn. OR YLIM EXCEEDED 5, /6H XL IM=E12.5,5X»6H YLlM=El2,5 // 19H WF 
PNOW GO TO BPLOT /// ) 
rFTURN 

20 xMAX= ALLOW (IT) 

TFIX=TI 
no 40 1 = 1 ,4 
1 1 = 1 

TF(YLTM *100. -ALLOW(T) ) 50 ,50 ,40 
40 CONTINUF 
60 TO 30 

50 yMAX rALLOW(TT) /10D, 
l CALL RSTFPM 

CALI GRIDGN < 1 23, 1 023, 24 ,Q24, 1 8 , 1 8*5, 5) 

CALL PLOT1 (1 ,1,ZFRO,XMAX,2ERO,YMAX»X,Y,nPLOT,1*1H/) 

Al 0NGX(U=0.0 
Al ONGY (1 ) =0 « 0 
no 60 1=1,6 

call labelx ( alongx(t) * i > 
call LAREI.Y ( aLONGY ( I ) , 1 ) 

A10NGX(I+I)= ALONGX(T) + .2* XMAX 
60 6lON6Y(I + n= ALONGY ( I ) + .2* YMAX 
CALL PRINT <200,975,1 2,0, 38 tXTTTLF) 

CALL PR INT(47, 200,0,12,38, YTITLF ) 

CALL PRINT(123, 1000, 12,0, 72, TTTLF) 
call DMPBIlF 

rFTURN 

END 


ooooo 
ooo to 

O0020 

00030 

00040 

00050 

00060 

00070 

00080 

00090 

OOIOO 

OOllO 

00120 

00130 

00140 

00150 

00160 

00170 

00180 

00190 

00200 

00210 

00220 

00230 

00240 

00250 

00260 

00270 

00280 

00290 

00300 

00310 

00320 

00330 

00340 

00350 

00360 

00370 

00380 

00390 

00400 

00410 

00420 

00430 

00440 

00450 

00460 

00470 

00480 
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SIBFTC BPLOT ROOOO 

FUBROl IT INF PPI OT (X*y*XLTM*YLIM* TITLE) POOlfl 

DIMENSION X ( 300 ) * Y ( 300) * YTITLF ( 1 0) » ALONGy(7) * YTITLE( io) 

niMENSION TITLE ( 12 ) P0030 

COMMON /ABC / ALLOW (7) * Al ONGX (7 > .NPLOT* ZFPO*XMAX, IFIX POOftO 

nATA (XTITLE(T) *I=t*lO)/3AH TIME (SEC.) / P0050 

DATA (YTITLE(T)*I=1»10)/3AH BONDLINF TEMPERATURE CR > / POOftO 

ALONGY (1 )=0,0 P0070 

no 10 1 = 1*7 POORO 

T I si P0090 

TF(YLTM -ALLOW( I > > 20*20,10 POlOO 

10 CONTINUF ROlin 

WRITE (6,1000) YLIM P0120 

1000 FORMAT (/// 37H PPLOT WILL NOT BE DONE BECAUSE YLlMs E12.S ////) P013O 

RFTUKN POIttO 

20 YMAX SALLOW (IT) R0150 

CALL RSTFRM P0160 

CALL OR I DON (123*1 023.24.924 *1B*1A*5*5 ) R0170 

CALL PLOT! (1,1* ZERO, XMAY *ZERO* YmAX*X * Y * NPLOT.1, 1H/ ) POlflO 

no 30 1=1*6 PD190 

CALL LABfcLX ( ALONGX(f) * 1 ) R0200 

CALL LABELY ( ALONGY ( T ) * 1 ) P0210 

30 ALONGY ( I +1) s ALONGY ( T ) + .2* YMAX P022ft 

CALL PRINT (200*975* 12* 0* 3**XTTTLE) P0230 

CALL PRIimT(47, 200*0*12. 3B.YTITLE) P02U0 

CALL PRINT (123*1000*12* 0 *72 »TTTl-F) R02R0 

CALL OMPhUF P0260 

PETURN P0270 

FND R02AO 
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SIBFTC CPLOT 

FUBROUTINF CPI OT (X, Y1 . Y?*XLIM,YI IM1 ,YLIm?,TITLF, Y) 

DIMENSION X(300) »Yl (300) ,Y2(300 ) .YTITlEMO) »YYC?000) ,XTITLE( 10) 
DTMFNSION TIT) E(1?)»Y(300)» ALONGY ( 7 ) 

DIMENSION CURVE(l) »VBH6CU) »HRH6<?) 

COMMON /ARC / ALLOW ( 7) * Ai.ONGX (71 .NPLOT * ZFPO, XMAX » IFIX 

oata (VRUR(I) ,i=i,4> / inn.o,Ro.n»20.o» sn.o / 

data (HRup(I) ,i=i, 7) / i.n»2.o,5.o#io.o»?o.n,50.o»ioo.n / 

DAT A ( XT ITLF ( I ) » 1=1 » 1 0) /.3RH TIME (RFC.) / 

oata (YTiTLF(T) ,I=l,lO)/38H DISTANCE (IN.) / 

DATA ONF/UH1 OF0 /.TWO/4H1A60 / 

DATA W0N/1H1 /,T00/1h? / 

C *** FOUR (4) CHARACTERS ARE ALLOWFD FOR CURVF(l) 

CURVE ( 1 )=ONE 
hFAC7R=HBUG( IFIX) 

FYMPOLSWON 

yPIR rAMAYl (VLlMIrYLTM? ) 

NCURVF =1 
no 1 T = 1 » NPLOT 

1 YV ( I )= Yl(I) 

no 7 T=1 »u 
Tl= I 

tF ( YBTG+ 100. -Al.LOW(T) )6.*,7 
7 CONTINUE 

WRITE (6,100n) YLIMI.YLIMR 

1000 FORMAT (/// 3<»H CPLOT WTU NOT PF DONF RFCAURP Y|_IMl=Fl?.S,tOH OP 
1 Yl IM2- FI?. 5 //// ) 

RFTURN 

6 yMAX =ALLOW (TD/100. 

\/FACTR=VBllG( IT ) 

CALL R5TFRM 

CALL ORIUON ( 1 23, 1 023, 24 .024 . 1 8 . 1 8 * 5# 5 ) 

.1=1 

70 no lo I=J, NPLOT 
TT = I 
IM =T-1 

tf( yy(d-yci) ) 20 »m#in 

10 CONTINUF 

N0PT=NPL0T-J+1 
LI =J + NOPT/? 

TVLOC=(YkiAX-YY(LL) )*1P.*\/FACTP +24. -4. 

THLOCs X (LL) *1 Q, /HFftCTP +123. -48. 

CALL PRINT (IHI OC.IVLOC, ft, 0»4, CURVE) 

CALL PLOT 1 ( 1 , 1 , ZERO , XMAX , 7ER0 , YMAX , X ( J ) , YY ( J ) , NOPT , 1 , SYMBOL ) 

TE (NCUPVF-i ) Po , 85 . PO 

85 no 86 1 = 1 »NP|.nT 

86 YY ( I ) =Y2 ( I ) 

CURVE U)=TWO 
PYMROLzTOO 

ncurvf = ? 

J=1 

fiO TO 70 
?0 NPT=IT-J 

LL=J + NPT /? 

IVL0C=(YMAX-YY(LL) )*1P.*VFACTR +24. -4. 

THLOCS X ( L L ) * 1 8 . /HFACTR +1?3. -48. 

CALI. PRINT ( IhI OC.IVLOC, 8,0 .4, CURVE) 


soooo 

sooio 

S0020 

S003P 

S004P 

soosn 

S0060 

S007P 

soobo 

songo 

sotoo 

sono 

soi?o 

SO! 30 
S0140 
SOISO 
SOI 60 
SOI 70 
SOI 80 
SOIPO 
S0200 
S0210 

so??n 

S023O 

S0240 

SO?SO 

S0260 

SO?70 

S0?80 

SO?PO 

S0300 

C0310 

S032O 

S033O 

S0340 

S03B0 

S0380 

S0370 

S0380 

C03PO 

S04fin 

«0410 

S04?0 

S0430 

S0440 

S04R0 

SOUftn 

S0470 

S0480 

S04PO 

SOSOO 

soPin 

SOS20 

S0530 

SOSUO 

S0S50 

sosfto 
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C at L PLO M <1 , 1 # ?FRO , XMAX . 7ER0 , YM AX , X ( J t » YY ( J ) ,NPT , 1 » 5YMP0L 1 50570 

DO 50 I J= II, NPLOT 50550 

JJ= I J 50590 

TF( YYt I«J)“ Y(TJ) )5n,b0*an 50500 

50 rONTIMUF 50510 

IF(NCllRVE-l) °0 * 85 * 9(1 50520 

40 J= JJ 50550 

GO TO 70 50540 

90 fil ONGYd )=0»0 50550 

no 100 1=1,5 «0550 

CALL LAREl.XUI ONGX(I) , 1) 50570 

CALL LAPEI Y(AI OM(iY(I),l) 50550 

100 A10NGY(I+1 )SA10NgY(I) + ,?*YMAX 50590 

CALL PRINT (200,975, 1 ?» 0, 55»XTTTLF) 50700 

CALL PRINT ( 47 ,200,0, 12*55 , YTITLF ) 50710 

CALL PRINT (125,1 000,1?, 0,72 »TTTLF) 50720 

CALL HMPHI IF 50750 

UPTURN 50740 

FA'O 50750 
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Fortran 

A 

AB 

ABLC 

ABLK 

B 

BB 

BL 

BLTEM 

BTEST 

C 

CB 

CHARC 

CHARK 

CP 

CPB 

CPC 

CFV 


APPENDIX C 
PROGRAM TERMINOLOGY 

Description 

"A" coefficient in matrix, single subscript 

"A" coefficient in matrix, double subscript 

specific heat of material at TABL 

thermal conductivity of material at TABL 

"B" coefficient in matrix, single subscript 

"B" coefficient in matrix, double subscript 

Total thickness of backup structure 

value of 1460 isotherm depth fran previous time step 

test to determine mode of heat transfer out of back surface 
of backup materials 

"C" coefficient in matrix^ single subscript 

"C" coefficient in matrix, double subscript 

specific heat of material at TCHAR 

thermal conductivity of material at. TCHAR 

specific heat of a node in ablation material 

specific heat of backup material node 

specific heat values in char specific heat table 

specific heat values in virgin specific heat table 

specific heat values in backup material specific heat 
tables 


CPX 
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Fortran 

D 

DB 

DELTT 

IMP 

DRH0 

DT 

UTS 

DX 

DXB 

D XV 

DXX 

EMBB 

EMC 

EMFB 

EMV 

EMX 

END 

EERl^ 

ERR2 

) 

ERR3 

ERR4^ 


Description 

"D" coefficient in matrix, single subscript 
"D" coefficient in matrix, double subscript 
time step in the time step table 

test used for dumping (IMP = 0 skip dump, IMP =1.0 
start dumping) 

local mass flow rate of ablation gas 
time step from the time step table in hours 
time step from time step table in seconds 
thickness of a node in the ablation material 
thickness of a node in a backup structure material 
variable ablation node thickness 

fixed ablation material node thi 

emissivity of back surface of each material in backup 
char material emissivity 

emissivity of front surface of each material in backup 
virgin material emissivity 

emissivity of front surface of ablation material 
code word for plot routine 

Control numbers for printing error statements when an input 
or calculational mistake is made 
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Fortran Description 

FBLjflW "blowing efficiency in reducing convective heating 

FC0NV factor to correct convective heating rate for various "body 
locations 

FENV emissivity — view factor product to cabin interior 

FRAD factor to correct radiative heating rate for various body 

locations 

FTEST test to determine mode of heat transfer into front surface 
of backup materials 

FV view factor for external environment 

G defined by Fortran statement 

GAPX gap width between backup materials 

H film coefficient between backup materials 

H300 enthalpy of air at 300° K 

HEAD any 72 alphanumeric characters used to identify problems 
being run - printed at top of first page of output 

HEADNG any 72 alphanumeric characters used to identify each input 
section 

HENV film coefficient to cabin environment 


HTX total enthalpy 

HV heat of degradation of virgin material 

HW wall enthalpy computed from enthalpy — temperature table 

HX enthalpy values in enthalpy table 

IEM test used to determine if front surface is virgin or char 

for vising proper emissivity 

IPRC variable print frequency in time-step table 

IFRCT present print control number 
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Fortran Description 

IE test to determine if node temperature is greater than TAEL 

IR1 test used in determining node density at TX1 temperature 

IR2 test used in determining node density at TX2 temperature 

NCASE number of problems to be run 

NCPB number of points in each backup material specific heat table 

NCPC number of points in char specific heat temperature table 

NCFV number of points in virgin specific heat temperature table 

HKC number of points in char thermal conductivity — temperature 

table 

HKPB number of points in each backup material thermal conduc- 

tivity table 

HKV number of points in virgin thermal conductivity temperature 

table 

HMB number of materials in backup structure 

HP number of node points in ablation material 

HPBS total number of node points in backup structure 

HPF total number of points in heat shield structure (HP + HPBS) 

HPL0T output plot control number 

HPM number of nodes per material in backup 

HHP number of points in enthalpy — temperature table 

HPTT number of points in time-step table 

HEEC number of points in surface recession — temperature or time 

table 

number of points in trajectory input table 


HTRAPT 
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Fortran 

RXA. 

NXB 

NXC \ 
KXD 

ME 

J 

QBL0CK 

QC0N 

QC0MX 

QHW 

QIN 

QL0SS 

Q0XID 

QRAD 

QRADX 

QUIT 

R 

RBI 

RB2 

RH0 

RH0BX 

RH0C 


Description 


dummy indexes for subroutine Save 


amount of convective heat blocked due to mass injection into 
boundary layer 

trajectory table convective heating rates 
cold •wall convective heat rate at present time step 
hot wall convective heat rate without blowing 
net heat flux into front surface 

boundary condition for heat transfer to cabin interior 

heating rate due to combustion 

trajectory table radiative heating rates 

radiative heat flux at present time step 

code word for plot routine 

thermal resistance due to conductivity between nodes in the 
ablation material 

thermal resistance due to conductivity between past and 
present node in backup material 

thermal resistance due to conductivity between present and 
forward node in backup material 

density of an ablation material node 

% 

density of individual materials in backup 
mature char material density 
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Fortran 

RH0V 

EH^fl 

RH0Y2 

S 

SD0T 
SAVEIT 
SAVEXX 
SAVX 
SAVY1 
SAVY2 
SAVY3 
SAVY4 
SAVY1X 
SAVY2X 
SAVY3X 
' SAVY4X 
SR 
T 

TABL 

TCHAD 

TCP 

TCPC 

TCFV 


Description 

virgin ablation material density- 

density of node at past time step 

density of node at present time step 

thermal capacity of a node in the ablation material 

surface recession rate 

depth of 1060 isotherm at any given time 

time corresponding to maximum depth of 1460 isotherm 

time corresponding to maximum depth of 1060 isotherm 

surface recession depth at maximum 1060 isotherm depth 

bondline temperature at maximum 1060 isotherm depth 

term that will contain maximum depth of 1060 isothexm 

depth of 1460 isotherm at maximum 1060 isotherm depth 

surface recession depth at maximum 1460 isotherm depth 

bondline temperature at maximum 1460 isotherm depth 

depth of 1060 isotherm at maximum 1460 isotherm depth 

term that will contain maximum depth of 1460 isothem 

surface recession values in surface recession table 

present time 

temperature at which ablation starts 
temperature at which ablation stops 

temperature values in backup material specific heat tables 
temperature values in char specific heat table 
temperature values in virgin specific heat table 
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Fortran 

TEMDI 

TEMPI 

TEW 

TEST2 

TDMP 

TIME 

TIET 

TITLE 

TKC 

TKV 

TL 

TLIM 

TREC 

TS 

TTABLE 

TTUL 

TUL 

TUL1 

TUL2 

TV 

T¥ 


Description 

arbitrary initial temperature distribution values 
constant initial temperature distribution value 
interior cabin temperature 

test to determine proper heat shield initial temperature 
distribution 

time to start dumping or printing information used in check- 
out of program (sets DMP = 1. 0) 

trajectory table time values 

starting time of problem 

control card used for reading in new data for successive 
problems 

temperature values in char thermal conductivity table 
temperature values in virgin thermal conductivity table 
total thickness of heat shield structure (VL + BL) 
time limit of problem 

surface temperature or time at which char removal is to 
start 

temperature or time values in surface recession table 
time values in time- step table 

equals TUL if VPT = 0 or equals TX2 if VPT = 1 — used in 
computing char properties 

maximum value of TX1 and TX2 

maximum TX1 values — used in computing gas ablation rate 
maximum TX2 values — used in computing gas ablation rate 
sink temperature of external environment 
temperature values in enthalpy table 
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Fortran 

TX1 

TX2 

TX2C 

TX2T 

TXK 

TX0 

TX 

VEL 

VELX 

VL 

VLI 

VLTEM 

VLV 

VPT 


WEKEEP 

XBM 

XC 

XI 

XIDET 

XK 


Description 

temperature of nodes at past time step 

temperature of nodes at present time step 

temperature at fixed locations in at la t ion material as de- 
fined by XC 

temporary storage of TX2 temperatures for computing thermal 
properties 

temperature values in backup material thermal conductivity 
tables 

initial temperature at front surface of heat shield for com- 
puting linear temperature gradient 

temperature distribution at forward time step 

trajectory table velocity values 

trajectory velocity at present tine step 

initial virgin material thickness 

initial ablation material thickness 

value of 1060 isotherm depth from previous time step 

variable ablation material thickness 

test to determine if properties are irreversible with 
temperature 

depth of 1460 isotherm at any time 
thickness of individual materials in backup 
fixed location of nodes in the ablation material 
node number 

any 72 alphanumeric characters to identify each material 

thermal conductivity values in backup material thermal 
conductivity table 
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Fortran 

XKB 

XKC 

XKV 

XL0ST 

XLSTI 

XLSTV 

XMDC 

XMDG 

XMD0 

XMDT 

XHP 

XNFM 

XPL0T 

XV 

YK 

YPL0T1 

YPL0T2 

YPL0T3 

YPL0TU 

ZZZ 


Description 

thermal conductivity of backup material node 

thermal conductivity in char thermal conductivity table 

thermal conductivity value in virgin thermal conductivity 
table 

amount of solid ablation material lost in a time step due 
to surface movement 

distance from original surface to present front surface lo- 
cation, inches 

distance from original surface to present front surface lo- 
cation, feet 

mass loss rate of char 

mass gas ablation rate due to pyrolysis of virgin material 

mass flux rate of oxygen to surface 

total ablation rate 

number of nodes in ablation material 

number of nodes per backup material 

time to be written on tape and plotted 

location of nodes in variable ablation material thickness 

thermal conductivity of a node in ablation material 

recession depth to be written on tape and plotted 

bondline temperature to be written on tape and plotted 

1060 isotherm depth to be written on tape and plotted 

1460 isotherm depth to be written on tape and plotted 

ratio to determine when the limiting value of heat blockage 
has been reached 
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APPENDIX D 
GENERAL FLOW CHART 



Continue 
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Go to initialization 
of program constants 
to start next problem 















TABLE I.- SAMPLE PROBLEM INPUT 


(a) Coding sheet 


. , M. .. J 


A 


JttPICKW C4MMIXH® KftUTCAR - T*$*\ C*V£ - 4/0/05 


AAiJ fill.;. 

40 QOOQ 


+ 0 © -+ 0 .© 
40 © 40 « I 
4 © © Hr © . V 


4* 00 
400 


t 

(00 
I 00 


I 


Al^aIi. , 40041 .0 




400 


. 44 ?] 


loaL, © ROte&XftO M W 


~J 1 I ■ 1, . ,i 1. , Inn 


1 4\000 ,;0 400 4V40©. o 
A©x.,05 1 400 4-0 . 75 

4 \.© 400 + 0.0 

•ii 2 2 9 

41 440 . O 400 40 . \ 2 
4141*0.. jO. +00 40. 43 
43^0. ou 4 oo 40 . oos 


L 4 A 00 i>'M 
ASOfl 
+300 
40 >0 


o 


+00 

400 

400 

400 


4 o.oC»(t 
+o » 0+3 
40 .43 
43.0 


40 o 

+00 

400 

J 

40 o 
+ 00 
400 
400 
+O 0 
400 

-04 


40 .© ... ,. 400 A 34 l. © 
4* Z9, 04+00i 4,1 aA_. 
40. \2 4*0 40.43 


L Z 
41 • o 
+ \ • o. 
+400-0 
+700. « 
4V©0O*O 

4» ( 00*0 
4000 * 0 


40.0 41.0 * A 

.400 . 42SOi. O 

400 40.07(0 , 40^ 


4i0.0. 4.0 « Oi 


400 , 


Aft*A 


+04 40 . tz 
4©4 +0»43 
400 40-0A5 


400 

400. 




AMI 


Am 


+0l.j43^. 400, 


ASAO. .A 4 ©O^ 4 . 0 .. 0055 - 401 . 


+00 +0 » 047 Z +00 *+30©.© .. 4-00 +.Q- I 04i+ 4 00 . , , j 
4(0 « 070 400 + \ \ 0© • 0 . 400 +0.070 400 , 


+ 0 O 

+00 

+00 


40 


.43 
. o 


r\ 


JXH- TRATeeT**\ - 

- . Z ' ■ 

40.0 400 435. o 


^*95 Rto/soo- 30^7 
4 Oo 40*0 • 400 +1*325 


+00 

-04 


+04 


( -A.-JI l_ 


> -U J ...~ 1 — J. . 


3 



. +00 ' ‘ 4* 5 . , O . 

* r 1 

■'liz 2zn 

+,S.«,o; i . . .4,©,© 

■ _ 3' i : 2, 

Mt #rxaa, 

4.3fo©«iQi l 4© © . . 4© • 9k S 
-±9Q © <»<* 

_ . V©« 40. ©<*9 

4© © k . 

. A ^Ai-tLQi- ]7_ 4.0 o 4© * \ 
+O.J 0 ! 7 L , +© o 40. © 


400 40-0 
400 


400 +2. *25 404 




.,ju4.:2,jla 




©,. \ xwcwe* 

40 0 4<(fc0.0 4,00 + 0* 0(*5 400 4540 1 0 400 4:0.0,45.5 4,00 

400 474 0.0 400 40.0411 400. 4140.0 40 0 40...0.42M 40.oj 

400 4 V040 . 0 4001 40 • 07 400 4\V4o.O 400 40. 0 7. . , 400! 

400 4 VVS0 ,© 40©J ,40.41 400 . ! 

40 o 40 . 3 .,... 400 4 o . 5 ^ 4 00 , , : 

400 40* o 400 40-0 400 




; ... MS0.4 i RK$teVr*K TV CftBTN 
L 454©, . 0 400 40 . o +00 

L ; T€H0«R<VTO*.C X* 

;_j4fiuj&L 400 455® .0 400 

) I I 

"7.4* 7 

, 40.0 400 40*0 400 

! 4Vi*l • 2 4.oo +44©o.o 4oo 
4\ V V 3 0 400 44000 .0 400 
1 4A400. o 4o© 44713. o +oo 
1 4 *700« o 400 4$ 2ft ft * © 4©o 


OWVIXR^HEIIT - 4«WV«0.0 
40.0 4 oo 4o*o 400 

Oo*5TR*T 
4530*0 400 


4342 . * 400 4 V4 00 . o 400 4449 . 7 4Q0. 4-Aft oo . © 4ftj _J 

4*m •» 4o© 43000 . 0 400 4943 -0 400 _43fefiO. 0 .40ft __ 

4\2oo. 0 4oo 44224.0 400 4A300.0 +00 44411*.. O 400 

4V500 O 400 4433<*. o 400 4l<*00.o 400 45V2TL.O iftft, 

4V$OO.o 400 45454.0 400 4*ft00. 0 400 +559W. o 400 



*1000.0 +oo +5713.0 +oo +3100.0 +oo +51SI .© +00 +2100 . o +oo + $9<»9*o +00 

I +2*00*0 +00 +<*011,© +00 + 2+00. 0 +00 +M3W.O +00 + 2500.0 +00 + <»29\.O+00 

+l(*oo.o +00 +439 S.o +O0 +H00.0 +00 +(0+91*0 +00 +U00.0 +00 ++S91 . o +00 

+2 + 0 0*0 +.00 +<•<.+$. 0 +00 +9000.0 +00 +(0905.0 -+00 +3VOO.O +00 +t+\$. O_ +_00 

+3100. 0 +00 +1050.0 +00 +3900.0 +00 +1V15.0 +00 +3+ 60 . 0 +00 +11S , Q-0-iOO 
+IS00 0 +00 +1+30.0 +00 +3+00*0 +00 +1V30.0 +00 +3100. 0 +00 +1AOO.O +ojb 

+ 33 00.0 +00 +1310.0 +00 +3300.0 +00 +9\ 10. 0 +00 ++Q0Q. o +00 +3360.0+66' 

+H\aO*o +00+950 0.0 +00 ++3 0 0.0 +00 +3100.0+00+^300.0 +00 +3136.0 +06" 

1 V*+'+aQ * o +00 +9000.0 +00 ++S00.0 +00, +91 50 . 0 +00 ++<000.0 +60 + 9110 *0+661 
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TABLE I. - SAMPLE PROBLEM INPUT 
(b) Fortran data card listing 


1 

lyf’ICAt Cl’AKPl AULAfOR - 7 FST C ASE - 4 / 6 /b 5 


+ r.nu o II 

+ llO 

+ r 

.u 

+un 

? 

1 

+ti.t) 

+ U n 

+ 1* 

.1 

+ nn 

100 


+ to f 'U«U 

+ u0 

+ ). 

.1 


1 UO 


+ 1 . u 

+o n 

+ i 

.0 

+un 




DONALO M. ClIRPY 


IYRTCAc CH#PKI-vO am. AT I ON matfrtal pkopfrtifs 


+ mb0.n 

+ll o 

+ 1 1 ibO . n 

+ 1)0 

+ D*U 

+oo 

+ 34.0 

+ 00 

+?o,o 

+on 

+ 0.00 

+ 00 

+ U .o5 

+ 0 0 

-M . 7«S 

• +on 


+oo 

+ 1.50 

+ 00 

+250,0 

+ 00 

+ 0.0 

+ 00 

f 1 .u 

41)0 

+ t*.0 

+ u 0 

+ D.1P 

+ 0O 

+ 0.43 

+ 00 

+0,070 

+ 00 

+ 0.43 

+00 


? 


2 ? 








+ 1 'u-, n . n 

4 ll 0 


+ U0 

+ 1.0 

+UO 

+ 0,1? 

+ 00 





+ I<‘b0 . 0 

4 (J 0 

+<* .43 

+ U0 

+ 1.0 

+0O 

+ 0.43 

+00 





4 360 . ii 

40 0 

+ i» » u* ^ 

+on 

+4^0.0 

+ 00 

+ 0 . U5f> 

+ 00 

+550.0 

+ 0 o 

+0.0555 

+00 

+»6u» u 

4 U 0 

+ i ! * U*o 

+ u o 

4/50,0 

+ 0 0 

+0.057? 

+ 00 

+850,0 

+on 

+0.0584 

+ 00 

+ 96 11 , 11 

41)0 

+ ti 9 U^H 

+ u0 

+ ioe>o,n 

+00 

+0.070 

+ 00 

+1150.0 

+oo 

+0.070 

+ 00 

+360.U 0 

40 0 

+ U .43 

+ 00 

+ iioo,n 

+uo 

+ 0,4.3 

+ 00 





4 b . U 

4 (| 0 

+ V.0 

-04 

+n00.0 

+00 

+9.0 

-04 





NO TkAj^CTuRY - 

0 =UR 

B1 H/ShC-SOFT 






C. 

4 (1 » U 

4-ijO 

4<-.6, 0 

+u 0 

+ 0.0 

+oo 

+?.9?5 

+ 04 





+t>n u • u 

+u n 

4MC..0 

+ U 0 

+ 0.0 

+oo 

+2.9?b 

+ 04 





1 


+ 0.1 

+ 00 









+ 3.0 

+ 11^ 











9 

? 











LSA( KLIP MATM>IAl 



0.1 

inches 

THTC* 




+ 3^ 0 • n 

40 0 

+ 0 . 0 6 5 

+ 0 o 

+450.0 

+00 

+0.055 

+oo 

+550.0 

+ 00 

+0.0555 

+ 00 

+ * U 9 U 

400 

+ 0 . 06b 

+ on 

+ /50.0 

+00 

+0.057? 

+nn 

+850.0 

+ 00 

+0.0584 

+ 00 

+ MMJ 9 0 

40 o 

4i ,069 

+ 1)0 

+ 1060.0 

+ 00 

+ 0.07 

+ 00 

+1160.0 

+ 00 

+ 0,07 

+ 00 

+3*0*00 

4||0 

4). .43 

+ 0 0 

+ 11 00. 0 

+00 

+0.43 

+ 0 0 





+J4 . 0 

40 0 

+ 0.1 

+ 1)0 

+0.9 

+ 00 

+ 0,9 

+oo 





+ U 9 U 

40 0 

+ H.0 

+ u0 

+ 0.0 

+ 00 

+ 0,0 

+ 00 





HP AT TRAiNSFFt TO TAP IN 

FmVIROMMFNT 

- HENV= 

0.0 





456(J . t) 

+ U ^ 

+ v . 0 * 

+ 1)0 

+ 0.0 

+ 1)0 

+ 0.0 

+ 00 





initial 

. iRMPt-RATgRE IS 

CONSTANT 







+u.o 

4(|0 

+h30.U 

+ 0 0 

+530.0 

+0 0 







U2 












+ u . 0 

4 0 0 

+ t>.0 

+ u0 

+302.9 

+00 

+1400.0 

+on 

+449,7 

+oo 

+1800.0 

+ 00 

4*>1 7.2 

4 0 0 

+2400.0 

+ 00 

+791.0 

+ 00 

+3000.0 

+oo 

+978,0 

+ 00 

+3500.0 

+ 00 

41113.0 

40 0 

+1+000,0 

+ 1)0 

+1POO.O 

+00 

+4?24,n 

+un 

+1300.0 

+ 00 

+4485.0 

+00 

4 1 u u o . n 

40 0 

4+1723,0 

4 0 0 

+1500.0 

+ 00 

+4Q35. 0 

+ 00 

+1500,0 

+on 

+5127.0 

+00 

41700.0 

400 

+b?99,n 

+ 00 

+ i«on.o 

+ 00 

+5U54 . n 

+oo 

+1 Q 00,0 

+00 

+5595.0 

+ 00 

4^00 0.0 

41 1 0 

+6728.0 

4 0 0 

+21 00.0 

+ 00 

+5551.0 

+00 

+2?00.0 

+ 00 

+5Q6P.0 

+ 00 

42^00 . 0 

41)0 

+ t-07B.O 

+ t)0 

+ 2.400,0 

+ 00 

+61 85.0 

+ 00 

+2500.0 

+on 

+6?91 .0 

+ 00 

4 261)0 . 0 

4 U 0 

4+395.0 

+ 00 

+2700.0 

+00 

+6497,0 

+on 

+2800.0 

+on 

+6597.0 

+on 

42°00.0 

40 0 

+e ? 699,n 

+ 00 

+3000. 0 

+00 

+bP05,n 

+ 00 

+3100.0 

+ 00 

+6918,0 

+ 00 
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f 32 n U .11 
+ ^ 0 (j . 0 

+ ^H0(| ,11 

+ 4 1 no . II 
1-U40u.il 


+ u n +7n c >u • u 
+ llO + 7u8|J,(| 
+ (l r, + 7'-»7u , (| 
+U O+ftKOQ . u 
+u r '+ Q iou , u 


+110+3300,0 
+ t|0 + 3oOO,U 
+uo+mno, 11 
+u0+420o,o 
-HiO+U'sno.n 


+ 1)0+7175.0 
+on+7b30.o 
+00+R1P0.0 

+on+s7no.o 
+ l)O+9l«M),0 


+ 00 + 34 00,0 
+00+3700,0 
+00+4000,0 
+00+4300. 0 

+no+uf,oo,n 


+0O+73S0.0 

+00+7800.0 

+00+8300.0 

+oo+P8so.n 

+00+P270.00 


+ 00 
+ 00 
+00 
+no 
+00 



TABLE II.- SAMPIE PROBLEM OUTPUT 


TYPICAL CHARRING ABLATOR - TEST CASJE - 4/6/65 _DONAi.D M. CURRY 


INPUT DATA. 


TIME L I M I T = 6 * OOOOE 02 INITIAL TIME=0. NPTT= 2 


TIME 

0. 

6. OOOOE 02 


TIME STEP 
1.0000E-01 
l • OOOOE-O J. 


PRINT CONTROL 
100 
100 


F dON V= _1 . 0QOQ 0E 00 FRAD= l.OOOOOE 00 


TYPICAL CHARRING ABLATION MATERIAL PROPERTIES 


TABL= 1.06000E 03 
FBLOW* 0. 

HV= 2.5'JOOOE 02 
CHARC= 4. 30000E-01 


TCHAR- 1.46000E 03 
EMV = 6.50000E-01 
VP T = 0. 

ABLK= 7.00000E-02 


TREC= 0. 

EMC* 7. 50000E-01 
F V- l.OOOOOE 00 
ABLC* 4.30000E-01 


RHOV* 3.40000E 01 
H300= l « 29060E 02 
TV= 0. 


RHOC* 2.00000E 01 
VL= 1.50000E 00 
CHARK* 1.20000E-01 


NP* . 31 NKC= ‘ 2 NCP C = 2 NKV= 9 NCPV= 2 NREC = 


TEMPERATURE 
3. 60000 E 02 
4. 60000E 02 _ 
5 • 60000 E 02 
6. 60000E 02 
7.60000E 02 
8.60000E 02 
9.60000E 02 


THERMAL 
CONDUCTIVITY 
6.5C000E-02 
6.5000QE-02 
6.55000E-02 
6.6C000E-02 
6.72000E-02 
6 • 84000E-02 
6.90000E-02 


VIRGIN MATERIAL.. 


TEMPERATURE 
3.60000E 02 
1 . 1 0000E 0 3 


SPECIFIC 

HEAT 

4.30000E-01 

4.30000E-Q1 


1.060001 0 3 7.0 OOOOE- 02 

l. 160006 03 7 • 00000 E- 02 


TEMPERATURE 
1 • 46000 E 03 
l.OOOOOE 04 

SURFACE RECESSION TABLE 


TIME SR - IN/SEC 
0. 9.00000E— 04 
6.000 00 E 02 9.000006-04 


CHAR MATERIAL 


THERMAL 
CONDUCTIVITY* 
1 . 20000 E- 01 
1 • 20000E-01 


TEMPE RATURE 


1.46000E 03 
l.OOOOOE 04 


SPECIFIC 

HEAT 

4.30000E-01 

4.30000E-01 


NO "TRAJECTORY - 0 = 95 "bTU/ SEC-SOFT 


NO. OF TRAJECTORY POINTS 


2 



TIME Q CONVECTIVE Q RADIATIVE VELOCITY 

0. 9.50000E 01 0. 2.92500E 04 

6 • OOOOOE 02 9.50000E 01 0. 2_._92500E 04 


PROPERTIES OF BACKUP '‘STRUCTURE 


NO. OF 

MATERIALS IN BACK-UP 

SHIELD* 1 




TOTAL NUMRER 

OF NODES IN BACK-UP SHIELD* 

3 



THICKNESS OF 

BACK-UP SHTELU* 

l.OOOOOE-Ol 

• . .... . 



BACKUP 

MATERIAL 

0.1 INCHES 

THICK 





THERMAL 



SPECIFIC 


TEMPERATURE 

CONDUCTIVITY 


TEMPERATURE 

HEAT 


3.60000E 

02 

6 . 50000E-02 


3.60000E 02 

4.30000E— 01 


4.60000E 

02 

6.5C000E-02 


1.10000E 03 

4.30000E-01 


5.60000E 

02 

6 • B5000E-02 





6. 600 00 E 

02 

6.60000E-02 





7.60000E 

02. . 

6 • 72000E-02 





8.6 0000 E 

02 

6 * 84000E-02 





9.60000E 

02 

6 • 90000E— 0? 





1 • O6OO0 t 

03 

7.00000E-02 





1. 16000 E 

03 

7 • COQ00E-02 






EMI SSI V I TY 


MATERIAL 

DENSITY 

THICKNESS 

FRONT 

BACK 

NOOES/MATERI, 

1 

3.4000E 01 

l.OOOOE-Ol 

9.0000E-01 

9.0000E-01 

3.0000E 00 


ADDITIONAL DATA FOR INDIVIDUAL MATERIALS" IN BACKUP STRUCTURE 


MATERIAL 

FILM COEFFICIENT 

GAP THICKNESS 

FTEST 

BTEST 

. _.L 

0. 

0. 

0. 

. o._ 


heat transfer to cabin environment - henv*o.o 

T E M^EKA T UiTlf* ~5 • 60000E~02 FTrM~OE F FIC I E NT* 0~. VIEW* FACTOR^ 6. 

INITIAL TEMPERATURE IS CONSTANT 


Q LOST* 0. 


TEMPERATURE DISTRIBUTION IN HEAT SHI ELD I S UNIFORM AND EQUAL TO 5. 3000 E 02 



t 


10U 

OUTPUT DATA* 


TIME- 9.90000E 00 QCONVECT I VE»_ 9, 5QOOOE 01 QRADI ATI VC* 0. VELOCITY* 2.92500E 04 

GAS ABLATION RATE* 0. CHAR ABLATION RATE* 5,40000E 00 TOTAL ABLATION RATE* 5.40000E 00 

RECESSION DEPTH* 9.00000E-03 QHOT WALL* B.8 9893 E 01 ; 

T EMPERATURE DIS TRIBUTION IN HEAT SHIELD AT THE END OF THE TIME STEP* T* l.OOOOOE 01 SECONDS 


TEMPERATURE DISTRIBUTION IN THE ABLATING MATERIAL 


4.25564E 

03 

2.70129E 

03 

1.19740E 

o? 

6.43320E 

02 

5.43814E 

02 

5 • 31490E 

02 

5.30126E 

02 

5 > 30007E 

02 

5 > 30000E 

02 

5.29999E 

02 

5.30000E 

02 

5.30000E 

02 

5.30000E 

02 

5.30000E 

02 

5.29999E 

02 

5.30000E 

02 

5 > 30000E 

02 

5* 30000E 

02 

5 >300006 

02 

5.30000E 

02 

5.30000E 

02 

5> 30000E 

02 

5.30000E 

02 

5 > 30000E 

02 

5.30000E 

02 

5.309C0E 

02 

5 > 300006 

02 

5.29999E 

02 

5> 30000E 

02 

5.30000E 

02 


5>29999E 02 


TEMPERATURE DISTRIBUTION IN THE BACK-UP STRUCTURE 


5.29999E 02 5.30000E 02 5.30000E 02 


TIME* 1.99000E 01_ OCUNVE CT I VE*_9 . 50000E 01. QRADJ AT I VE* _0 . VELOCITY* 2*92500E 04 

GAS ABLATION RATE* 0* CHAR ABLATION RATE* 5.40000E 00 TOTAL ABLATION RATE* 5.40000E 00 

RECESSION DEPTH*' Jl« 600 OOE-P2 QHOT WALL* 8.89148E 01 ; 

TEMPERATURE DISTRIBUTION JN HEAJ SHIJLD AT THE E NDOF THE_ T I ME STEP, T* 2. 00000 E 01 SECONDS_ 


TEMPERATURE D I ST R I BUT ION I N TH E~ A R L A TING MATERIAL 


A.26944E 03 _3>15689£ 03 

5.37192E 02 5.31229E 02 

5.29999E 0? _ 5.29999E 02 

5 * 29999E 02 5.29999E 02 

. 5 _»2_9 9 9 9 E — 02 5. 2 999 ? ,E.,.g2 

5.29999E 02 


1.77750E 03 
5.30183E 02 
5 # 29999E _02 
5.29999E 02 
5.29999E 02 


1 *02157E 03 
5* 3002 3 E 02 
5.29999E 02 
5 • 2 9999 E 02 
5> 29999E 02 


6.89430E 02 
5.30002E 02 
5.29999E 02 
5.29999E 02 
5.29999E 02 


5 • 66334E 02 
5* 30000E 02 
5.29999E 02 
5.29999E 02 
5 >29999£ 02 


TEMPERATURE DISTRIBUTION IN THE BACK-UP STRUCTURE 


5.29999E 02 5.29999E 02 5»29999E 02 


TIME* 2 >990006 0 1 QCON VECTIVE* 9>50000E 01 QRADI ATI VE= 0. _ __ VELOCITY* 2 *92 500 E 04 

GAS ABLATION RATE* 4>04922£ 00 CHAR ABLATION RATE* 5.40000E 00 TOTAL ABLATION RATE* 9.44922E 00 
RECESSION DEPTH* 2.70000E-02 QHOT WALL* 8.87279E_01 _ 
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Figure X. - Radiation temperature approximation error. 
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Figure 2.- Schematic diagram of charring ablator thermal protection system. 
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Figure 3.- Thermogravimetric data for typical charring ablation material 
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Figure - Charring Material property variation used as input 
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figure 5- - Comparison of temperature histories for nanahlating steel slat (pure conduction). 
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Figure 6.- Comparison of temperature histories for moving boundary model, 


no 
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Figure 8. - Plot program surface recession curve from typical charring ablator test case, 
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Figure 9. - Plot program bondline temperature curve from typical charring ablator test case 






